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SUMMARY 


An experimental program to develop grack growth and fracture 
toughness data under mixed mode conditions was undertaken. In a 
unique mixed— mode machine, room-temperature tests were made of 21 
flawed aluminum plates in the shape of a 91x91-cm (36x36~in.) mod- 
ified Maltese cross. Six were 2219-T87 1.55-cm (0.61- in.) thick; 
four were 2219-T87 0.64-cm (0.35-in.) thick; and eleven were 7075- 
T7351 1.27-cm (0.50-in.) thick. Ten specimens were tested stati- 
cally to fracture, five were tested cyclically then statically, 
and six were only tested cyclically, under loads ranging from pure 
tension on the flaw tip to pure shear. All flaws were oriented 
45° to the loading directions. Analysis of test results was di- 
rected to computing critical strain energy release rates, G , and 

strain energy release rate, AG, versus crack growth rate, Aa/AN. 
Stress intensities were computed from coefficients derived from 
finite-element analyses of straight, Z-shaped, and branched cracks. 
The Maltese cross shape permitted any combination of K^. and K 

with a 45° flaw orientation. Stress intensity coefficients were 
nearly constant for a range of crack lengths.- The area of strain 
uniformity under biaxial load in the center of an unflawed speci- 
men extended for a diameter of about 25.4 cm (10 in.). 

During cyclic testing, it was impossible to maintain a high 
proportion of shear-mode deformation on the crack tips. The cracks 
either branched or turned, resulting in crack fronts with little 
shear. Under static loading to fracture, straight 45° cracks re- 
mained straight only when shear stress intensity exceeded normal 
stress intensity, otherwise they turned. 

Mixed-mode crack-growth rate data, analyzed in terms of G, 
compared reasonably well with published single-mode data. There 
was little difference in the effectiveness of analyses used to 
analyze crack growth data. Straight and branched cracks were 
analyzed using Holston factors, and Z-shaped cracks using Iida 
factors, yet there was not an abrupt change in computed quanti- 
ties when cracks changed configuration. Measured crack displace- 
ment during cyclic tests agreed well with straight and branched 
crack analyses over the range of applicability. 

In fracture tests, values of G for pure shear were approxi- 
mately 50% higher than G for pure normal opening, and there was 
a large reduction in K^. resulting from application of'Kj. T . How- 
ever, net section stresses were well into the inelastic range when 

fracture occurred with high shear. Values of G were affected by 

cr J 

the relation between flaw growth direction and rolling direction 
for the 7075-T7351 alloy but apparently not for 2219-T87. 



INTRODUCTION 


Technical Background 

The development of linear elastic fracture mechanics has pro- 
vided concepts and data to designers that may be used to formulate 
criteria for the structural design of elements and the selection 
of materials. These criteria are based on characterization of 
flaw growth in the element by parameters that describe the stress 
field near the ends of the flaw. The specification of this stress 
field for an arbitrary loading condition requires determination of 
the three stress intensity parameters K^., K ^ and K^^.. These 

parameters provide a convenient means of relating the stress field 
conditions to the crack, loading conditions. Loading conditions 
such that displacements are normal to the crack plane (opening 
mode) are represented by the parameter K^. Loading conditions 

producing shear displacements parallel to the crack plane and in 
the plane of the specimen are represented by K^, and loading 

conditions producing shear parallel to the crack plane and per- 
pendicular to the specimen are presented by K . 

ft 

Most fracture mechanics work, both analytical and experimental 
has focused on the opening mode of deformation. Tensile speci 

mens were loaded normal to the flaw; analyses were developed for 
different configurations with loadings that produced only the open- 
ing mode. However most structures experience at least biaxial load 
ing. For example a spherical tank has equal biaxial stresses and 
a cylindrical pressure vessel has a 2 to 1 stress state. Further- 
more the flaw may be oriented at any angle relative to the prin- 
cipal stresses since the flaws are uncontrollable. Thus the open- 
ing mode is a special case, with the combined or mixed mode being 
the general case. 

In the past the mixed-mode problem has been ignored in most 
designs. The flaw has been assumed to be oriented perpendicular 
to the largest principal stress. Then allowable flaw sizes have 
been established based on the critical value of and flaw growth 

has been calculated considering only K . This approach has been 

unavoidable due to the lack of understanding of the mixed-mode 
problem and the corresponding absence of appropriate data. How- 
ever, there is no reason to believe that it is necessarily con- . 
servative . 



Only limited data are available for definition of the effects 
produced by the introduction .of shear loading that is concurrent 
with the opening mode load. Information developed during a recent 
NASA program, NAS3-14346 (ref. 1) , has revealed that a significant 
reduction in the tensile load-carrying ability of a cracked 4340 
steel member may accompany the application of mode II loading con- 
ditions. Limited data on some 5000- and 7000-series aluminum 
alloys obtained from the literature show much less degradation 
than xtfas found for the 4340 steel alloy. The degradation effect 
of combined modes I and II mixed-mode loading appears- to he ma- 
terial-dependent for static fracture conditions and virtually 
undefined for cyclic conditions. Crack* growth data — both static 
and cyclic — for materials of structural interest in aerospace 
hardware are needed to evaluate the consequences of high shear 
loading conditions. The continued lack of data concerning "the 
mixed-mode fracture behavior will continue to block use of frac- 
ture control measures in the design of aerospace hardware and 
may lead to serious structural deficiencies. 


Ob j ective 

The objective of this program was to develop crack growth 
and fracture toughness data under combined loading "mixed-mode" 
conditions. To accomplish this objective aluminum plates with 
through flaws have been cyclically and statically tested under 
biaxial loading and the resulting data analyzed*.- The purpose of 
this work was also to establish the applicability of various 
theories of flaw growth and fracture to mixed-mode conditions. 

To this end the test results have been compared with predictions 
gained from the methods -of analyses. 


Scope 

The objectives were pursued through a program of room tem- 
perature testing in a unique machine of 21 flawed , aluminum plates 
in’ the- shape of a 91.4x91.4-cm (36x36-in.) modified Maltese cross. 
All flaws were through the thickness, 5- to* 10-cm (2- to 4-in.) long 
and oriented 45 deg to the directions of load application. Three 
combinations of alloy and thickness were employed — six specimens 
were 2219-T87 1.55-cm (0.61-in.) thick, four specimens were .2219- 
T87 0.64-cm (0.25-in.) thick, and 11 specimens were 7075-T7351 
1.27-cm (0.50-in.) thick. ' 



Ten specimens were tested statically to failure, requiring 
one to three minutes to fracture* The mixed-mode combination of 
loads ranged from pure tension on the flaw tip to pure shear 

( K Il) # specimens were tested cyclically until substantial 

crack growth was realized and then tested statically to fracture. 
The six remaining specimens x*ere tested only cyclically. Although 
the cyclic tests were also intended to span the full range from 
pure K_£ to pure K^, curved crack .trajectories and branching phe- 
nomena resulted' in essentially removing shear from the crack tips. 

The analysis of the test results consisted, of (1) computing 

critical stress intensities. K-,_ and , and critical strain 

* Ic lie 

energy release rates, G , for the static tests based on plots of 

the crack growth resistance, R, and (2) plotting strain energy 
release rate, G, versus crack growth rate, and compliance versus 
crack length for all the cyclic tests. The test results are 
compared with predictions derived from the theories of Griffith 
as modified by 'Irwin (extended), Erdogan and Sih (maximum stress), 
and Sih (strain energy density) . Variables taken into account 
in the analyses and comparisons were alloy, thickness, final roll- 
ing direction, and the elffect of the xatio of Test re- 

sults are also compared with data, from other investigations. 



STATE-OF-THE-ART REVIEW 


Combined-Mode Fracture Data 

The study was limited to determining the effects on static 
fracture and fatigue of combined modes I and II loading conditions. 
For static fracture conditions, the effects of combined inodes I and 
II loading are commonly described in terms of an interaction dia- 
gram such as the one shown’ in figure 1. Figure 1 plots data from 
tests by a number of investigators using several materials in at- 
tempts to define modes I and II interaction effects ♦ The data 
range from indicating little or no interaction effect for static 
fracture to indicating a very significant interaction effect. 
Additional data of a similar nature for 7075-T6 and 7075-T651 indi- 
cate almost no interaction effect in one case (ref. 2) and a fairly 
significant amount in another (ref. 3), although not as strong as 
the linear interaction seen by Shah (ref. 1) and Liu (ref, 4). 

Ranganath and Goolsby (ref. 3) reported a significant effect 
on interaction of the material thickness . This implies that plas- 
ticity effects could play an important role in the interaction proc- 
ess . In fact, .plastic zone sizes are quite different in modes I 
and II. The associated energies to produce fracture could then be 
quite different and the interaction diagram might depend on the way 
in which the modes I and XI plasticity fields interact for struc- 
tural metals. 

It appears that the K^. and calculations at fracture in most 

if not all of the referenced data assumed that the crack did not 
undergo out-of-plane growth before reaching instability and the 
K^, values plotted did not account for the out-of-plane crack 

growth that normally occurs in mixed modes I and II fracture situ- 
ations. In this study, subcritical growth was experienced and at- 
tempts are made to present interaction data in which the effects of 
subcritical growth and growth direction were taken into account. 


Theories of Fracture 

i 

Table I summarizes the various theories of fracture, some of 
which are intended to predict fracture under conditions of mixed- 
mode loading. Included in the table is information on the physi- 
cal basis for the theory, the predicted ratio between mode I and 
mode II critical stress intensities, the predicted crack propaga- 
tion direction, and the predicted interaction diagram. 



Source 


Material 


Shah (ref. 1) 

4340 Steel 

Liu (ref. 4) 

7075-T651 A 

Liu (ref. 4) 

2024-T3A 

Wilson (ref. 18) 

7178-T651 A 

Pook (ref. 19) 

A 

A 

DTD 5050 A 

A 

- 

“V s- Lower Bound from 
X Wilson & Pook Data 


Figure X. - Interaction of mixed-mode loading on plates 







TABLE I. - SUMMARY OF FRACTURE THEORIES 


Identification 

Basis 

K_ t /K t 
lie Ic 

Fracture direction 

Mode I and Mode II Interaction 

Griffith 

Strain energy release rate 
G = G 

cr 

Not applicable 

Assumed colinear growth 

Not applicable G a G — — «1 

cr G 

cr 

Irwin/Griffith 

Related Kj. to G 
1-v 2 

G a — Kj = G (plane .strain) 

Not applicable 

Assumed colinear growth 

’ , K I 

Not applicable Kj. a K Ic y~ “ 1 
. * Ic 

Extended Irwin/Griffith 

Added shearing mode energies 
G “ +■ Kj-^J = G cr (plane strain) - i 

1.6 

Valid only for colinear growth 

t/Z v Z 

K 2 + K 2 = K 2 — L + — “ 1 

. i 11 Ic- 

Maximum stress theory 
Erdogan, Sih (ref, *1 ) 

Maximum tangential stress intensity 
CT ee v/KF= cos I ( K i cos2 1 ' 1 k ii sln 6 ) = K ic 

0.86 6 

Normal to maximum o aA 
* 0u 

° c - slirl (j sin *s) - 

0 / K I V 6 c 

cos COS 2 r - . 


reduces to Irwin/Griffith for “ 0 

' 

/K 

• . -! 5 sln M = 1 

■ ' 



h 

/ 

Strain energy density 
Sih (ref, 8) 

Minimum strain energy density 

u - f \ - k/vt 

S “ s or “ a U k l 2 + 2a 12 k r k 2 + a 22 k 2 2 
Reduces to Irwin/Griffith for « 0 


Direction of minimum S 
2(l-v) sin (0 c - * E ) 

- 2 sin [ 2 ( 0 C - ^ 

- sin 2e D 0 
' ' c 

/Ik 

+ 2 a 12 sin cos <$ e + 

+ &2Z cos2 = 1 

• 



-J 1 
* 11 

|3-4v - cos 6 J fl + cos ej 
a H “ 4(l-2v) 





3,, - 2 8l “ 0 [ cos 9 c-( 1 “ 2v >] 
4(l-2v) 

• 



' 

a„ - 4(1 - u) ( 1 - cos 6 c) 

4(l-2v) 





(1+cos $ c ) ^3 cos 0^ - 
4(l-2v) 




The first theory identified is the Griffith theory of fracture. 
This fracture theory is based on the concept that unstable crack 
propagation will occur when the strain energy released during an 
increment of crack growth is bigger than the energy necessary to 
drive the crack through the material. In this theory the energy 
required to drive the crack was considered to be surface energy 
associated with the new crack surface. Cracks are assumed to grow 
in the original plane of the crack, which is oriented normal to 
the applied tensile stress. In its originally stated form, this 
theory does not apply to situations involving mixed-mode fracture. 

The second theory is identified as the Irwin/Grif f ith theory. 
The physical basis for this theory is the same as the physical 
basis for the Griffith theory, i.e., the strain energy release 
rate is. considered to be the driving energy source for the fracture 
process. However, in this theory the energy required to drive the 
crack through the material is- considered to be dominated by the 
energy of plasticity required to form a small plastic zone at the 
crack tip as the crack propagates through the material. This 
energy is thought ;to be very much larger than the surface energy 
of the Griffith theory. For conditions of highly constrained 
plasticity, a material property G cr is considered to characterize 

the fracture process. This theory also introduces the notion of 
stress intensity factor and develops a relationship between the 
stress intensity factor and the strain energy release rate G by 
considering the amount of energy released as a crack in an elastic 
field grows by an incremental amount in the original plane of the 
crack. The fracture criterion is then recast by considering a 
critical value of stress intensity at failure known as fracture 
toughness. Uiis theory also does not directly relate to conditions 
of mixed-mode fracture. 

.The third theory referred to is the extended Irwin/Griffith 
theory. This theory attempts to deal with problems of mixed-mode 
fracture- on the same physical .basis as used in the first two 
theories. In this case, however, the total strain energy release 
rate for a mixed-mode situation is computed by adding the strain 
energy release rates associated with the modes of fracture to be 
considered. It is important to note that adding the strain energy 
release rates due to the .separate modes in this way is valid only 
as long as the crack is assumed to undergo incremental growth in 
the original plane of the crack. Fracture is still considered to 
occur when the strain energy release rate reaches a critical value. 
Under this assumption, the ratio of mode II fracture toughness to 
mode I fracture toughness is found to be unity. The predicted in- 
teraction is given by the equation shown for that case in table I, 
which is the- equation of a circle in the first quadrant of a Kj, 

space. 


9 



The extended Irwin/ Griffith theory can be made applicable to 
situations of mode X, mode II loading where crack growth occurs in 
a plane other than the original plane of the crack. To do this it 
is necessary to know in advance*, or be able to predict by an appro- 
priate hypothesis, the angle at which crack growth* will occur. The 
stress intensity factors appropriate for the new crack growth di- 
rection must then be' calculated. In this manner the calculation of 
G through adding the -strain energy release rates due to the active 
modes is compatible with the assumptions of the theory. This ap- 
proach has been discussed by Cotterell (ref. 5) , Hussain, et al. 

(ref. 6), and others. 

The fourth theory to appear in table I is identified as the 
maximum stress theory discussed by Erdogan and Sih (ref. 7). The 
physical basis for this. theory is the hypothesis that fracture 
occurs along a radial line emanating from the crack tip on which 
the tangential stress component is a maximum. It is further hy- 
pothesized that fracture occurs when the stress intensity factor 
for tangential stress along that line reaches a critical value, 
which is taken to be the fracture toughness of the material nor- 
mally associated with a pure mode I test. This theory reduces to 
the Irwin/Griffith theory for a pure mode I situation. A ratio of 
mode II to mode I fracture toughness of 0.866 is predicted by this 
theory, and this value does not depend on the mechanical properties 
of the material. The predicted fracture direction is given by the 
equation shown in table I and depends on the ratio of applied mode 
II stress intensity to mode I stress intensity. The equation for 
the predicted interaction is also given in table I and, when plot- 
ted in a space, has a shape similar to a parabola in the 

first quadrant. 

The fifth theory shown is identified as the strain energy den- 
sity theory due to Sih (ref. 8). The theoretical basis of this 
theory relates to calculation of the density of strain energy in 
the region near the tip of a crack loaded in mixed mode. The strain 
energy density is found to increase in' a singular fashion as the 
crack tip is approached along a radial line. The intensity of this 
singularity is identified as S and is assumed to take on a critical 
value at fracture. Further, the direction of crack propagation is 
hypothesized to be predicted by the direction in which S takes on a 
minimum value. As the other theories described, this theory reduces 
to the Irwin/ Griffith predictions for a pure mode I loading situa- 
tion. The ratio of mode XI fracture toughness to mode I fracture 
toughness predicted by this theory is shown in table I and is seen 
to depend on the elastic properties of the material. The equa-j 
tion that predicts the direction of crack growth is also given in 
the table and depends on the ratio of to K^. The equation for 

mode I, mode II interaction given in the table is found to again 
appear much like a parabola when plotted in the first quadrant of 
Kf, space. 



Figure 2 presents a plot of the angles of crack propagation 
predicted by the various theories and shows that, in the range of 
interest, there is little difference between the predictions of 
crack angle. Figure 3 is a plot of the interaction diagrams that 
are predicted by the various theories. Again, it is clear that 
there is very little difference between the theories in predicting 
mode 1, mode XX interaction. 


* Mixed-Mode Fatigue 

Very little information is available concerning fatigue crack 
propagation under mixed-mode loading conditions . This is because 
cracks, when subjected to mixed-mode conditions, have been found to 
rotate their direction of propagation immediately on cyclic load 
application in such a way that they propagate in primarily a mode 

I condition. Xida and Kobayashi (ref. 9) considered mixed-mode 
crack propagation in 7075-T6 aluminum. The mixed-mode condition 
was obtained by placing a' straight crack in a tension- loaded strip 
at an angle to the loading direction. Various degrees of mixed- 
mode loading were obtained by placing the crack in the plate at 
varying angles . It was found that immediately on onset of cyclic 
loading the crack began to turn in a way that would tend to make 
the crack plane perpendicular to the loading direction. This pro- 
duced a crack configuration- similar to that shown in Figure 4 taken 
from reference 9. 

To properly analyze the test results, Iida and Kobayashi per- 
formed a finite element stress analysis of cracks with geometries 
similar to those shown in figure 4. Figure 5 shows the results of 
that analysis. The analysis clearly shows that as the crack grows 
along its path, out of the plane of the original crack, the stress 
intensity situation remains mixed mode for only a short interval 
after which the loading condition is primarily mode I. In connec- 
tion with mixed-mode crack growth rates, Iida and Kobayashi con- 
cluded that the presence of mode II accelerates crack growth rate 
by 10 to 20% over that expected by a mode I fatigue test. 

Roberts and Kibler (ref. 10) have also studied problems of mode 

II fatigue crack propagation. In their experiments, cracks were 
loaded in combined mode I and mode II where the mode II component 
of the load was cycled but the mode I component of the load was 
held steady. Under these circumstances it was possible to produce 
crack growth in the original plane of the crack. However, in some 
cases branching of the crack xras observed. It was generally found 
that an increase in the steady, mode I load for a given cyclic 
a mp litude of mode II variation increased the crack growth rate. 






a 



a 


Initial Crack Orientation . 

^Figure 4.' - Crack extension pattern of an initially slanted crack. 

Recently,. Yokobori, et al. 9 (ref. 11) have conducted a series 
of torsional fatigue experiments on cracks in thin-walled hollow 
cylinders'. The^material used for these studies- was 5083P-0 alumi- 
num alloy. The experiments were conducted in a way that subjected 
the crack to alternating pure shear stress with no normal stress 
present. Under these loading conditions, the crack was found to 
branch at each of the crack tips as diagrammed in figure 6. This 
configuration implies that the branches of the growing cracks are 
subjected primarily to the mode I stress intensity factor during 
the growth portions of loading for each of the branches. 

From the scant amount of information available on mixed-mode 
fatigue crack propagation, it is evident that little is understood 
about this phenomenon. Because the mode of crack growth in most of 
the experiments conducted to date has been primarily mode I, it has 
been possible to interpret available data in terms of the tradition- 
al mode I fatigue laws. There has also been sufficient evidence 
that the effect of cyclic mode II in the presence of mode I loading 
is "to accelerate the crack growth rate, possibly by significant 
amounts. This effect makes further study of this type of phenomenon 
important. " 
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Normalized Stress Intensity Factor, M $1 = K^V^or M S j j * Kj 



0 0.1 


Normalized Projected Crack Length 

Figure 5. - Normalized stress intensity factors for a tension plate with 
initially slanted crack. 
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Biaxial Specimen Stress Analysis 


The specimen used in this program was analyzed by IIo Is ton 
(ref. 12) using finite elements with mixed-mode crack tip singu- 
larity elements. The specimen geometry, including the crack, is 
shown in figure 7 and the normalized mode I and mode XI stress 
intensity factors are- given in figure 8. The flaw half-length 
scale used in figure 8 is shown in proportion to the specimen size 
in figure 7. In this program the tab width was 50.8 cm (20 in.). 


t t t t t t t 



until 
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Figure 7. - Specimen geometry analyzed by Holston. 
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EXPERIMENTAL PROGRAM 


Specimens 

Materials . - Two aluminum alloys were investigated in this 
program — 2219 in thicknesses of 1.55 and 0.64 cm (0.61 and 0.25 
in.) and 7075 in a thickness of 1.27 cm (0.5 in.). The material 
was received in the form of plates 1.219 m (4 ft) by 3.658 m (12 
ft). Each plate was cut into four specimen blanks 0.925 m (3 ft) 
square and four coupon blanks 0.3084 m (1 ft) by 0.925 m (3 ft). 
The thinner 2219 material was received in the T-87 condition. 

The thicker 2219 material was received in the T-37 condition, was 
cut, then aged at 177°C (350°F) for 18 hours to T-87. The 7075 
was received as T-651, was cut, then aged at 163°C (325°F) for 
28 hours to T-7351. Room temperature tension tests were made on 
30.5-cm (12-in.) long full-thickness flat dogbone coupons, with 
a gage section 5-cm (2-in.) long by 1.27-cm (^~in.) wide, cut 
from the coupon blanks, two coupons for each direction from each 
alloy/thickness combination. There was one longitudinal strain 
gage on each coupon. The coupon test results are given in table 
II. 


Preparation . - Fifteen 2. 54 -cm (1-in.) diameter holes were 
.drilled along each edge of the specimen blank for attaching grip 
plates, and four 2.54-cm (1-in.) diameter holes were drilled to 
serve as the roots of notches later put in with a band saw. The 
plan form of the specimens is shown in figure 9.* They were full 
plate thickness. A 0.635-cm (%-in.) diameter hole was drilled 
through the center of the plate to permit the insertion of a saber 
saw blade. A through flaw 5 to 10-cm (2 to 4-in.) long was then 
cut with a saber saw 45 deg to the loading axes of the specimen. 
The flaw width was 1.2 mm (3/64 in.). For the last 1% mm ,(0.05 
in.) at each end of the saw cut, a special blade was used, made 
from two pieces of band saw blade welded together and then ground 
to a V. A steel razor blade was pulled through the flaw at each 
end to establish a sharp notch for starting a crack on the first 
two specimens. This practice was later abandoned because it did 
not improve the precracking behavior. After flaw cutting, both 
surfaces of the specimen were polished in the area around the 
flaw using an orbital sander with 200 grit, then 600 grit, emery 
cloth followed by buffing with a cloth wheel and jewelers' rouge. 
The specimen was then ready for flaw sharpening by precracking. 


*Strain distributions for this shape are described in Appendix C. 
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PLATE MATERIAL 


Young’s modulus 

Ultimate strength 

Elongation 
at failure, 

% 

^m 2 x 10 6 

psi x 10 6 

M/cm 2 x 10 3 

psi x 10 3 

7.4 

10.8 

48.2 

69.9 

12.0 

7.4 

10.7 

48.3 

70.1 

10.0 

7.4 

10.8 

48.3 

70.0 

11.0 

7.6 

11.0 

48.1 

69.7 

9.0 

7.4 

10.8 

47.8 

69.4 

9,0 

7.5 

10.9 

48.0 

69.6 j 

9,0 

7.2 

10.4 

51.4 

74.5 

12.2 

7.2 

10.4 

• 51.4 

74.6 

12.2 

7.2 

10.4 

51.4 

74.6 

12.2 

7.3* 

10.6 

51.7 

75.0 

11.2 

7.2 

10.5 

51.6 

74.8 

11.0 

7.3 

10.6 

51.6 

v 74.9 

11.1 

7.2 

10.5 

46,9 


10.0 

Zi3 . . 

10.6 

47.1 

68.3 

10.0 

7.2 * 

10.5 

47.0 

68.2 

10,0 

7.2 ' 

10.5 

47.1 

68.3 

11.0 

7.2 

10.4 

47.2 

68.5 

11.0 

7.2 

10.4 

47.2 

i 

68.4 



11.0 
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Precracking . - Precracking to sharpen the flaws was per- 
formed by loading the flaxes hydraulically. A fixture was de- 
signed and built to accomplish this in a manner similar to that 
reported by Abbatiello and Derby (ref. 13). The precracker 
consisted of two 25xl5-cm (10x6-in.) 4130 steel plates 5-cm 
(2-in.) thick as backup plates for two 2.4-mm (0.093-in.) thick 
0-ring retainer plates. The backup plates were held in place by 
61-cm (2-ft) long clamping beams, one on either side of the test 
specimen (fig. 10 and 11). Hydraulic pressure was applied in one 
of two manners. Originally .the hydraulic power supply of a 5-kip 
MTS testing machine was used at 7 cycles/ second with the actuator 
bypassed. This limited the pressure to 2400 N/cm 2 (3500 psi) . 
Later, hydraulic pressure was supplied by a hydraulic ram that 
was placed in the testing space in the testing machine and at- 
tached to the loading plates. In this way the ram became a 
pressure intensifier and pressures up to the design limit of the 
ram and hoses could be obtained (6900 N/cm 2 or 10 000 psi) al- 
though these levels were never needed. Development trials to 
determine pressure levels and number of cycles to obtain 1% to 2% 
mm (50 to 100 mils) of crack growth were run to act as a guide 
for precracking the test specimens. These trials and the pre- 
cracking parameters for the specimens, including the final flaw 
lengths, are presented in Appendix B. , 


Test Equipment 

The equipment especially fabricated or gathered into systems 
for this program included the specimen grips, the testing machine, 
a COD gage holder for biaxial measurements, a compression fixture, 
and recording and data analyzing systems. 

Grips. - Flat plate grips were used to transfer the load from 
the testing machine clevises to the specimens. Each set of grips 
•consisted of eight plates, one for each side of each loading tab, 
with holes at one end to match the loading holes in the specimens 
and a large hole at the other end to accommodate the loading pin. 
Two sets of grips were used. . A steel set was fabricated first to 
permit transferring the full capability of the testing machine 
into the specimens. When it was determined that the required 
load values were much less, an aluminum set was made to reduce 
Weight and .increase the efficiency of specimen installation and 
removal. The aluminum grip plates are shown attached to a specimen 
in figure 12 . 




Figure 10. - Precracker O-ring retainer plates. 



Figure 11. - Precracker in place on specimen. 
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Figure 12. - Aluminum grip plates attached to specimen. 

Testing machine . - A special biaxial testing machine, the 
mixed-mode testing machine (MMTM) , was built for this program. It 
comprises a 1 million pound MTS universal testing machine, the 
vertical axis, to which has been added a specially supported 
horizontal axis, both under computer control. Both axes are 
operated off the same 0.53-mVmin (140-gal/min) hydraulic power 
supply. Figure 13 is a view of the MMTM. Figure 14 is a sche- 
matic diagram of the horizontal axis. It fits between the four 
columns of the MTS machine and is located so the test specimen 
can be mounted between the two axes. The weight of the hori- 
zontal axis is reacted by a support system designed to minimize 
the "crosstalk 11 between axes. Functionally the horizontal axis 
is identical to the vertical axis. The frame reacts the load be- 
tween the hydraulic actuators and the passive end whiffletree. 

A schematic diagram of the horizontal axis support system is shown 
in figure 15. By ultimately resting on compressed gas, the hori- 
zontal axis can move up and down with the horizontal centerline of 
the specimen without loading the specimen with its dead weight. 

The rams behave essentially as very soft springs. They are pinned 
at top and bottom to eliminate any longitudinal restraint of axial 
deformation of the horizontal axis frame. There is some loading 
of the specimen during cyclic operation of the vertical axis due 
to the inertia of the horizontal axis. This is a small amount and 
is discussed in Appendix C, Mixed-Mode System Checkout Summary 
Report, where the entire checkout procedure results are given 
because proper functioning of this system was so vital to the valid- 
ity of the program test results. 
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Figure 13. - Mixed-mode testing machine. 
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Figure 15. - Schematic diagram of horizontal axis support system. 


The horizontal axis has its own load cell circuitry and power 
supply, balancing circuit, output amplifier, fail-safe circuits, 
and independent control of the servo valves. Thus it can be oper- 
ated independently of the vertical axis, either manually or under 
computer control, or locked in with the vertical axis under com- 
puter control. The horizontal axis load cell and its circuitry 
was calibrated by comparing its output with that of the vertical 
axis, using a calibrated "load link” as an intermediary. First 
the MTS machine was calibrated in the standard fashion by the 
Denver Division Metrology Laboratory. Then an extensively strain- 
gaged, dogbone-shaped aluminum load link was inserted in the 
vertical axis. The gage output was read on a portable strain 
indicator at various load levels. Next the link was placed in 
the horizontal axis. Now acting as a ‘‘secondary standard,** it 
was used to calibrate the output of the horizontal axis load cell. 
Excitation voltage and amplifier settings were established to give 
a 10-V output for full-range loadings, thus matching the MTS 
machine. Shunt resistors, a permanent part of the load cell cir- 
cuitry, were used to check the stability of the output over long 
periods of time. 


A PDP 11/05 Digital Equipment Corporation computer with 8k 
memory was used to control the mixed-mode machine for all test- 
ing, both cyclic and static. A two-channel arbitrary program 
exercised real-time control throughout a test run, making the 
computer a sophisticated function generator. Forty values of 
load versus time can be put on each axis before the program 
repeats, with the signal between values being either a ramp or 
a haversine. On this project the ramp function was always used 
and no more than 10 values were ever needed to define the cyclic 
load shape. Frequency of loading was not limited in the com- 
puter but rather by the flow rates of the hydraulic fluid, mass 
of the servo valves, elasticity of the structural components of 
the machine, etc, and depended on the load levels and faithful- 
ness of response desired. The cyclic rate was kept below 30 
cycles/minute because of the soft response of the horizontal 
axis, generally testing at a rate of 15 to 20 cycles/minute. 

Biaxial displacement gage holder . - The analysis in Appendix 
A predicted both opening and shear deformations of the flaws in the 
specimens. A special device was designed and built to measure 
these deformations independently at the center of the flaw. It 
was a split cylinder that fits inside a ^>-in. diameter hole at the 
flaw center. The cylinder halves incorporated knife edges and bear- 
ing wheels and were slotted to permit opening displacement and com- 
pletely reversed transverse displacement without interference. 

Motion was detected by two double cantilever strain-gaged displace- 
ment measuring devices. One acted against a knife edge and bearing 
wheel located at the middepth of the specimen. It held the split 
cylinder in place while the adhesive dried and recorded the normal 
crack opening displacement (COD) . The other acted against a knife 
edge and a wheel that protruded from the specimen and recorded crack 
shear displacement (CSD) . Figure 16 presents a view of the dis- 
assembled cylinder and the two measuring devices. Figure 17 shows 
an in-plane view of the holder and gages mounted on a specimen. 

The normal mode opening measurement side is shown in figure 18 
and the shear mode side in figure 19. 

A new double cantilever gage was built for this program to 
match the performance of the Martin Marietta-built device shown 
in figure 19, except that the cantilever leaves were narrower than 
those shown in the figure to facilitate insertion into the mid- 
depth of the split cylinder. The biaxial gage assembly was oper- 
ationally verified by mounting it into a plate cut in two to 
permit free motion of the parts in a testing machine or on the 
table of a traveling microscope. Motion was imposed in one prin- 
cipal direction while readings were taken in the other to deter- 
mine whether there was "crosstalk . 11 Motion of 0.76 mm (0.030 
in.) in the "normal opening" direction resulted in motion in the 
"shear" direction of about 0.0025 mm (0.0001 in.) or 1 part in 
300. 
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Figure 16. - Biaxial displacement gage holder 
disassembled. 



Figure 17. - Biaxial displacement gage assembly in 
place, in-plane view. 
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Figure 18. - Biaxial displacement gage assembly, 
normal opening side. 



Figure 19. - Biaxial displacement gage assembly, 
shear mode side. 


ORIGINAL PAGE® 

OF POOR QUALITi 



Compression fixture . - Some of the tests required compression 
to be applied on one axis of loading. A fixture was designed and 
built to permit the application of compressive loading on the 
vertical axis while supporting the specimen and the grip plates 
against lateral buckling. Figures 20 and 21 are two views of a 
specimen ready for testing with the compression plates in place. 
Figure 22 is a schematic of the system. The system consisted of 
two compression plates, one on either side of the specimen, clamped 
to the grip plates at one end and free to slide in a slot at the 
other. The compression plate on one side was clamped at the op- 
posite end from the compression plate on the other side. Each com- 
pression plate had two 10-cm (4-in.) wide 2.5-cm (1-in.) thick bars 
welded to it in the vicinity of the notch roots of the specimen. 
These bars rested against the specimen and provided lateral sup- 
port while permitting in-plane deformation of the specimen. Thus 
all axial load introduced into the grip plates by the testing 
machine was transmitted to the specimen except for the slight load 
carried by friction in the slots or between the specimen and the 
lateral support bars. These surfaces were lubricated before each 
test. The compression plate had a hole in the center to permit 
movies to be taken of crack growth and to provide access for the 
crack displacement measuring gages. Buckling occurred in only 
one test (SN-12) out of the eight in which the plates were used. 

Monitoring and recording equipment . - Data were monitored and 
recorded by several redundant sets of equipment. Three Dana model 
5600 digital voltmeters were used to continuously monitor the hori- 
zontal axis load cell and two load cells in the vertical load train, 
the MTS load cell, and the strain gages on the bottom pull rod. 
Instrumenting the bottom pull rod permitted it to be used as 
another load cell to provide a check on the performance of the 
horizontal axis support rams in keeping the horizontal axis weight 
off the test specimen. These digital voltmeters were also used to 
monitor the specimen installation and removal procedure and the 
setting of preloads when required. A 25x43-cm (10xl7-in.) X-Y 
plotter was used as a diagnostic tool during the cyclic tests and 
to record loads during the static tests. It was generally used 
for recording the vertical axis load versus the horizontal axis 
load. During the cyclic tests, it was thus possible to tell if 
the loads bore the proper time phase relation to each other and 
whether any adjustments had to be made in peak values by fine- 
tuning span settings on the control consoles. During the static 
tests, the plotter provided a permanent record of the fracture 
loads. A teletype keyboard was used for inputting the two-channel 
arbitrary programs into the computer and for periodically record- 
ing the load cell feedback signals. It also made a permanent 
typewritten record of the number of cycles of loading and the 
elapsed time for each block of cycles. A cycle counter on the 
machine control console provided a running total of cycles of 
loading. 



Figure 20. - Compression fixture installed on specimen, console side. 



Figure 21. - Compression fixture installed on specimen, opposite side. 
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Figure 22. - Schematic diagram of compression fixture attached to specimen 



Both crack displacement gages, both load cell outputs, both 
loading piston displacements, and time of day were intermittently 
recorded on a Honeywell model 1612 visicorder during the cyclic 
tests and continuously during a static test. The displacement 
gages were excited by an independent power supply. All seven 
signals were amplified and conditioned in a Martin Marietta-built 
cart prior to being fed to the visicorder. The cart also con- 
tained calibration switches for all the channels. Two oscillo- 
scopes were used to monitor whatever signals happened to be of 
interest at any time. 

Subcritical crack growth and fracture were recorded with a 
D. B. Milliken model 5-2 camera. This camera has a variable fram- 
ing speed of 1 to 500 frames/second with an accuracy of ±1%. It 
accommodates a 122-m (400-ft) roll of 16mm film. Film speed can 
be changed while the camera is running and two channels of time 
data (specific event blips and/or clocks) can be put on the film 
edges. It was run at 100 f rames/second . Figure 23 is an enlarge- 
ment of a frame from the motion pictures taken at 100 frames/second 
during the static test of specimen SN-18. The numbers are from 
light-emitting diodes attached to Fluke model 8000A digital volt- 
meters driven by the output of the load cells on each axis. The 
circles are scribed on the specimen at 0.254-cm (0.1-in.) inter- 
vals to aid in crack tip location. Crack growth versus load was 
established by projecting the film with an L and W photo-optical 
data analyzer, model 224-A, and measuring the projected crack 
length, suitably scaled, at each increment of load displayed on 
the digital voltmeter readouts in the picture. 

A running historical log kept of each test recorded date, 
time of day, notes on any anomalies, nominal peak load values and 
cycle rates, number of cycles of loading for each load block and 
cumulative number of cycles, scribe mark number (see following 
section), and the symbol used to tie together the teletype record, 
visicorder record, and X-Y plot at any given time (usually an 
alphabetical letter). In addition, film speeds, static loading 
rates, observations on specimen behavior, and changes in record 
scale factors were recorded. 


Test Procedures and Parameters 

Specimen installation . - The first step in preparing a pre- 
cracked specimen for test was to mount the grip plates. A 2.4x 
2.4-m (8x8-ft) plywood assembly table with appropriately spaced 
standoffs was built to align the specimen and grip plates during 
assembly. Grip plate bolts were tightened with an impact wrench 
to a 135 to 200 N-m (100 to 150 ft-lb) torque. If the compres- 
sion fixture was required, it was installed along with the grip 
plates. Next the 0.64-cm (%-in.) center hole, used as a starter 
for flaw cutting, was enlarged to 1.28 cm (h in.) by drilling and 
reaming and the crack displacement gage holder was bonded in place. 
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Figure 23. - Frame from high-speed motion picture of specimen SN-18 
during static testing. 


A scribe mark was made across the ends of the sharpened flaw and 
labeled 0. The horizontal axis was rolled out from between the 
vertical axis columns and the specimen with grip plates attached 
was lowered by crane between the horizontal axis side frames, 
proper note being taken of the final rolling direction and flaw 
direction. Shear pins were slipped through the grip loading holes 
and clevises on the horizontal axis load rods with appropriate 
shim plates to keep the specimen midplane on the load train center- 
line. The horizontal axis with specimen in place was then rolled 
back between the vertical axis columns. The support rams were at- 
tached and the horizontal axis was lifted about 3.8 cm Oh. in.) 
off its rails. Then the vertical axis shear pins were slipped into 
place, care being taken not to load the specimen by adjusting the 
support, rams while monitoring both the upper and lower vertical 
axis load cells. 

Amplifier gains on the crack displacement channels were ad- 
justed before each test. Before they were inserted into the bi- 
axial gage holder, the gages were exercised in a supermicrometer 
and the gains set to provide a convenient excursion of the traces 
on the visicorder. Calibrations on all the channels were checked 
with built-in shunt resistors. The computer was programmed accord- 
ing to a prepared test plan and the resulting computer output 
signals were verified on the X-Y plotter. Range set tins, load 
limits, error limits, plotter and visicorder zero settings, etc 
were all reviewed. These settings were reviewed again and ad- 
justed as needed after the hydraulic power supply was switched 
to the test pressure level and the loading pistons were adjusted 
to the zero load position. - 

Load types and directions . - The objective of pursuing a test 
program that covered the range from pure K^. to pure conditions 

on the flaw required several different loading schemes. It is 
theoretically possible to obtain the full range of ^ rom 

zero to infinity on a 45-deg flaw by using the cyclic loading 
schemes shown in figure 24 provided the flaw remains at 45 deg to 
the loading axis while it grows. All the schemes shown (except 
E) were designed to reverse the direction of shear, or on 

each cycle. It was believed that reversing the shear on each cycle 
would cause the flaws to grow in a stair-step fashion on the micro- 
scopic scale, and that they would therefore appear straight on a 
gross scale. For the static tests, the scheme designation applies 
only to the first ramp of the load. No cyclic tests were made 
using scheme D (see Table III for the load type applied to each 
specimen). Because the crack branched or turned, the intended 
ratios were often not obtained. 
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TABLE III 


TEST PARAMETERS 


Spec- 

imen 

Alloy a 

Thick- 
ness, 
cm (in.) 

Final 
rolling 
directxon 
in test 

machine'* 

Grip 

mate- 

rial 0 

Intended 

nominal 

k iAx 

Test Type 

load 

scheme 

Initial Flaw Lengths 

Cyclic 

Static 

2a i? cm 
(in.) 

2a^, cm 
(in.) 

2a*, cm 
(xn.) 

bA/6, 
cm (in.) 

SN-1 

2219 

1.55 

H ! 

S 

Pure K_ 

Yes 

No 

A 

6.541 



8.64 



(0.61) 



1 




(2.575) 



(3.40) 

SN-2 

2219 

1.55 

V 

s 

Pure K t . 

Yes 

Yes 

A 

5.385 

12.55 

10.29 

8.15 



(0.61) 



I 




(2.120) 

(4.94) 

(4.05) 

(3.21) 

SN-3 

2219 

1.55 

V 

A 

2.4 

Yes 

No 

B 

5.469 






(0.61) 







(2.153) 




SN-4 

2219 

1.55 

H 

S 

1.0 

Yes 

Yes 

C 

5.423 


9.09 

5.423 



(0.61) 







(2.135) 


(3.58) 

(2.135) 

SN-5 

2219 

.0.64 

II 

A 

2.4 

Ye3 

Yes 

B 

5.283 


7.82 

5.283 



(0.25) 







(2.080) 


(3.08) 

(2.080) 

SN-6 

2219 

0.64 

H 

S 

Pure K t 

Yes 

Yes 

A 

4.928 

12.27 





(0.25) 



1 


A 


(1.940) 

(4.83) 



SN-7 

7075 

1.27 

V 

s 

Pure K_ 

Yes 

No d 

A 

5.166 



9.09 



(0.50) 



J. 




(2.034) 



(3.58) 

SN-8 

7075 

1.27 

H 

A 

2.4 

Yes 

No 

B 

6.375 






(0.50) 





a 


(2.510) 




SN-9 

7075 

1.27 

V 

A 

2.4 

Yes 

No d 

B 

5.474 



5.474 



(0.50) 







(2.155) 



(2; 155) 

SN-I4 

7075 

1.27 

V 

A 

2.4 

Yes 

Yes 

B 

5.029 


7.52 

5.84 



(0.50) 







(1.980) 


(2.95) 

(2.30) 

SN-27 

7075 

1.27 

H 

A 

Pure K tt 

Yes 

No 

E 

7.701 


5.44 

7.701 



(0.50) 



II 




(3.032) 


(2.14) 

(3.032) 

SN-U 

2219 

0.64 

H 

A 

0.4 

No 

Yes 

D 

10.60 

10.62 





(0.25) 



i 




(4.175) 

(4.18) 



SN-12 

2219 

0.64 

H 

A 

Pure K II 1 

No 

Yes 

E 

10.73 

10.73 

7.59 

10.732 



(0.25) 






(4.225) 

(4.23) 

(2.99) 

(4.225) 

SN-15 

7075 

1.27 

II 

A 

Pure 1^ 

No 

Yes 

A 

10.26 

10.26 





(0.50) 






(4.040) 

(4.04) 



SN-1 6 

7075 

1.27 

H ' 

A 

1.0 

No 

Yes 

C 

10.92 


7.72 

10.922 



(0.50) 







(4.300) 


(3.04) 

(4.300) 

SN-17 

2219 

1.55 

H 

A 

0.4 

No 

Yes 

D 

10.49 

10.49 





(0.61) 







(4.130) 

(4.13) 



SN-1 8 

2219 

1.55 

H 

A 

Pure K II 

No 

Yes 

E 

10.34 

10.34 





(0.61) 






(4.070) 

_(4.07) 



SN-23 

7075 

1.27 

H 

A 

1.0 

No 

Yes 

C 

10.57 


7.47 

10.566 



(0.50 







(4.160) 


(2.94) 

(4.160) 

SN-24 

7075 

1.27 

H 

A 

Pure K_„ 

No 

Yes 

E 

.10.52 

10.52 





(0.50) 



II 




(4.140) 

(4.14) 



SN-25 

7075 

1.27 

V 

A 

Pure 

No 

Yes 

E 

10.17 

10.17 





(0.50) 






(4.005) 

(4.01) 



SN-26 

7075 

1.27 

H 

A 

0.4 

No 

Yes 

D 

10.72 

10.72 





(0.50) 







(4.22) 

(4.22) 




a 2219 , temper T-87; 7075, temper T-7351. 

- horizontal; V - vertical. 

C S - steel; A - aluminum. 

^Last cyclic load taken as a "static" test. 
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Direction of final rolling , - Care was taken to note the 
relation between the final rolling direction and the orientation 
of the specimen in the testing machine/ This relation was de- 
liberately and frequently changed in the early stages of the 
program to eliminate any bias that might be due to always putting 
the rolling direction in the same axis. Later in the program it 
was changed to study the effect of rolling direction on the path 
of crack growth , fracture direction, and fracture loads. Final 
rolling directions are listed in table III. 

Flaw size . - There are several values of !, initial n flaw size 
associated with each specimen. The first is the length of the 
flaw after precracking, designated 2a^. The next is the arc 

length of curved flaws after cyclic testing but before static 
testing for ^specimens subjected to both types of tests and used 

for Holston r s analysis ^designated 2a|?j. There is also the pro- 

jection of this length in the direction perpendicular to the 
final fracture direction used for Iida T s analysis and designated 
I I 

2a^. The designation 2a^ is also used for the projected length 

of branched cracks analyzed by the Ilda method for fracture tough- 
ness. Finally,, there is the equivalent linearized slant length 
for the Ii'da analysis used to establish the relation between the 
crack length and the plate width, designated 6. Figure 25 

contains diagrams illustrating these various concepts of flaw 
length, as well as an illustration of the length designations 
used in the crack growth analyses. The values are given in Table 
III along with the other test parameters. 

Cyclic test procedure . - The first step in performing the 
cyclic tests was to program the computer. It was assumed that 
both loading axes would behave as though infinitely stiff and 
follow the computer output exactly. A few cycles of load, usually 
less than *10, would then be applied and the loads recorded on the 
visicorder at a high rate of paper speed (10 cm/s). The difference 
in the time each axis reached its peak load could then be read 
easily to 0.01 seconds. Changes were made in the computer pro- 
gram for the vertical axis to delay it so the two axes would 
reach the peak loads simultaneously. This delay ranged from 
0.03 to 0.07 second depending on the magnitude and rate of load- 
ing. The X-Y plotter was used to record vertical versus hori- 
zontal load before and after the program change. Figure 26 is 
an example of the plots obtained in this manner for SN-1 at the 
565-kN (127-kip) load level. 




(I Ida Analysis) 






After 

Precracking 


(Holston Analysis) 


For Curved Cracks 


For Branched Cracks 


Starting Flaw Lengths for Static Tests 
, 4 2a (iida) 



Note : T = top, p 2a. 

B = bottom, \ 4a . R ✓ i 

V = vertical, \/ X 

FI = horizontal. 'e 

Holston analysis a Vertical crack 2a = 2a. + 4a WT + 4a x/D , 

I V ! Vd 

horizontal crack 2a « 2a. + 4a UT + 4a UD . 

I HI Hb 

Iida analysis = Vertical crack 2a = 1.414a. + 4a VT + 4a VB> 
horizontal crack 2a = 1.414a. + 4a UT + 4a UD . 

I HI Hb 
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Horizontal Load, P u 

H 



Figure 26. - Example of computer program change to overcome 
horizontal axis lag. 

Cyclic load was applied in accordance with two different 
schedules. In some tests it was applied in 7 to 10 groups of 
three load blocks each. Each load block consisted of a number 
of cycles necessary to produce about 0.13 cm (0.050 in.) of flaw 
growth. Each of the three blocks in a group was run at approxi- 
mately the same maximum load level although the stress intensity 
increased as the flaw length increased. Each group of three 
blocks was run at increasingly higher load levels, the first 
group at a stress intensity of about 25% of critical, with the 
maximum load being increased by 10 to 20% for each subsequent 
group. In other tests the load level was increased 5 to 10% 
after each load block. The ratio of minimum to maximum loads 
(R value) were always kept below 0.1 and generally below 0.05. 
Cyclic loading was generally continued until the crack growth 
rate reached 0.127 to 0.38 mm/ cycle (5000 to 15000 micro-in./ 
cycle). In some cases it x-zas continued to fracture. The cycle 
rate x*as generally 15 to 20 cycles per minute. In plates of 
finite extent, and' are different functions of crack lengthy 

even though the crack stays straight (see Appendix A, fig.' A-ll). 
Therefore to maintain a constant ratio of as the crack 

grew it was necessary to slightly adjust the ratio of the load on 
one axis to the load on the other. Figure 27 contains curves of 
the ratios of loads as a function of half crack length for a 




45-deg crack required to keep the K y/ K ratio constant ♦ This 

chart was used as a guide even though the cracks often turned or 
branched, making the analysis not entirely valid. At the end of 
each load block, the locations of the tips of the crack were 
scribed on the specimen and noted in a data book. A section of 
oscillograph record was also taken to permit determination of the 
slope of a crack opening versus load plot, ttfhich is a measure of 
the compliance. Load from the vertical axis versus load from the 
horizontal axis x*as recorded on an X-Y recorder at the beginning 
of each load block to check the load ratios. The computer print- 
out and a cycle counter on the vertical axis console kept track 
of the number of cycles of loading at each load level. If the 
specimen was to be tested statically it was left in the machine. 
Otherwise it was removed and those that were still in one piece 
were either pulled apart or cut with a saw to expose the crack 
surface and to facilitate measurement of the distance between 
scribe marks. 



Half Crack Length, a, in. 

Figure 27. - Test load ratio chart. 
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Static test procedure . - The first step in performing the 
static tests again was to program the computer. There was no 
problem of horizontal axis lag because of the slow rate of load- 
ing. The rate used was 448 to 667 kN (100 to 150 kips) per minute 
for the axis with the higher load. The ratio of loads depended . 
on the desired ratio. The program was carried well above 

the expected maximum load. The high-speed movie camera was mounted 
and focused, and was turned on after some loading had taken place 
but before subcritical crack, growth started. It was run at 100 
frames/second to the end of the test. Oscillograph records were 
taken continuously throughout the test as was a record of loads 
on the X-Y plotter. When fracture occurred the error signals 
were generally great enough to shut off the machine, although in 
some instances failure was not violent and the machine was stopped 
manually. In some instances the crack displacement gages were 
removed before^ fracture to prevent damage to them. After the- 
test was over, the specimen was removed by reversing thie procedure 
for specimen installation, except care 'bad to be taken in unpinning 
the broken specimens as they were now generally in two distinct 
parts . 


Experimental Observations and Primary Data 

Because the type of analysis applied to the response of 
each specimen depended on that response, it is necessary to de- 
scribe how the cracks grew and under what conditions. Therefore, 
a short description of the behavior of each specimen under test 
is given in Appendix D, which also contains tables of cyclic load- 
ing schedules and resultant crack lengths for all cyclic tests, and 
tables of loads and crack lengths taken from the high-speed movies 
for the static tests. Load and cycle data are very briefly sum- 
marized in table IV. Appendix E contains photographs of all speci- 
mens after testing. Characteristics of crack appearance and other 
pertinent notes about the specimens are summarized in table V. 

Several general observations should be made. Great care 
was taken to align the loading axes of the testing machine and to 
locate the midplane of the specimens on the plane defined by these 
axes through the use of shims and spacers. Nevertheless there 
Ttfas evidence of specimen bending in almost all the tests. Some- 
times the specimens appeared to move ,out of plane at the center 
as much as 0*6 cm ft in.), while at other times no such movement 
was discernible although there was other evidence of bending. In 
the tests of specimens under pure conditions there should 

have been no response from the displacement gage measuring shear 
deformations (CSD) . Yet in these tests (SN-1, 2, 6 and 7) the 
CSD did have some output from the very beginning of the tests. 



TABLE IV. - SUMMARY OF LOAD AND CYCLE DATA 


Specimen 

Maximum 

Total 

Load at start 

Load at 

static 


cyclic 

cycles 

of subcritical 

fracture, kN 


load, kN 


crack growth. 

(kip) 


(kip) 


kN (kip) 






P 

P 

P„ 

P TT 




V 

H 

V 

H 

SN-l 

575.6 

52 263 










(129.4) 







- 

SN-2 

591.6 

59 689 

594.3 

604.1 

770.0 

779.3 

SN-3 

(133.0) 

489.3 

22 654 

(133.6) 

(135.8) 

(173.1) 

(175.2) 






(110.0) 






. 


SN-4 

569.4 

31 970 

745.5 

0 

849.6 

44.5 


(128.0) 


(167.6) 

(0) 

(191.0) 

(10.0) 

SN-5 

264.7 

43 105 

323.4 

137.9 

545.4 

230.4 


(59.5) 


(72.7) 

(31.0) 

(122.6) 

(51.8) 

SN-6 

243.8 

55 005 

297.1 

306.5 

508.0 

518.7 

SN-7 

(54.8) 

511.5 

40 455 

(66.8) 

(68.9) 

(114.2) 

(116.6) 






(115.0) 


SN-8 

23 673 





513.8 






(115.5) 


SN-9 

39 306 





502.6 






(113.0) 








SN-14 

313.6 

124 310 

1023.0 

0 

1468.0 

0 

SNr-27 

(70.5) 

435.9 

26 326 

(230.0) 

(0) 

(330.0) 

(0) 






(98.0) 








SN-11 

SN-12 

— : 


-153.0 
(-34.4) ' 

-315.8 

371.0 
(83.4) 

314.0 

-256.2 . 
(-57.6) 

-411. 0 b 

614.7 

(138.2) 

404. 8 b 



SN-15 

— 


(-71.0) 

667.2 

(70.6) 

667.2 

C-92.4) b 

965.3 

(91.0) b 

967.5 



SN-16 

(150.0) 

0 

(150.0) 

794.0 

(217.0) 

0 

(217.5) 

1108.0 






(0) 

(178.5) 

(0) 

(249.0) 

SN-17 




a 

-613.9 

1486.0 



CL 

(a) 

-676.1 

cl 

(a) 

680.6 

SN-18 

(-138.0) 

-1103.0 

(334.0) 

1130.0 



SN-23 

(-152.0) 

689.5 

(153.0) 

0 

(-248.0) 

751.7 

(254.0) 

0 



SN-24 

(155.0) 

-547.1 

(0) 

547.1 

(169.0) 

-1034.0 

(0) 

1025.0 



SN-25 

(-123.0) 

-578.3 

(123.0) 

578.3 

(-232.5) 

-1008.0 

(230.5) 

1001.0 



SN-2 6 



(-130.0) 

a 

(130.0) 

a 

(-226.5) 

-498.2 

(225.0) 
1281. 0° 



CL 

CL 




(a) 

(a) 

(-112.0) 

(288. 0) C 


Undetermined. ^Buckled. c Loading tab failure. 
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.TABLE V. - SUMMARY OF CRACK CHARACTERISTICS 


Specimen 

High-speed 

movies 

Crack growth characteristics 

Remarks 

1 SN-l 

No 

Curved to horizontal. 

No records of fracture. 

SN-2 

Yes 

Curved to vertical. 

— 

SN-3 

No 

Branched s horizontal pre- 
dominant. 

Cyclic test stopped- 

SN-4 

Yes 

Branched, horizontal pre- 
dominant at fracture. 

Tension normal to rolling 
direction in static test. 

SN-5 

Yes 

Branched, horizontal pre- 
dominant at fracture. 

High tension normal to rol- 
ling direction in static test. 

SN-6 

Yes 

Straight, then curved as 
shear lip formed. 

— 

SN-7 

No 

Curved to vertical. 

Cycled to fracture. 

SN-8 

No 

Branched, horizontal pre- 
dominant . 

Loading tab fracture during 
cycling. 

SiI-9 

No 

Branched, vertical predom- 
inant. 

Cycled to fracture. 

SN-14 

Yes 

Branched, horizontal pre- 
dominant at fracture. 

Shear not reversed on each 
cycle. 

SN-27 ' 

No 

Directly vertical with no 
curvature 

Cycled to fracture. Shear 
not reversed. 

SN-11 

Yes 

Straight into corners, no 
shear lip . 

Tension parallel to rolling 
direction. 

SN-l 2 

Yes 

Sharply curved to vertical. 

Buckled under vertical com- 
pression. 

SN-15 

Yes 

Gradual curve to horizontal. ; 

— 

SN-16 

Yes 

Vertical, then horizontal 
into corners . 

Tension parallel to rolling 
direction. 

SN-l 7 

Yes 

Straight into corners, no 
shear lip. 

— 

;SN-18‘ 

Yes 

Straight into corners, no 
shear lip. 

— 

SN-2 3 

Yes 

Horizontal. Compare with 
SN-16. 

Tension normal to rolling 
direction. 

SN-2 4 

Yes 

Straight into corners, no 
shear lip. 

Tension parallel to rolling 
direction. 

SN-2 5 

Yes 

i 

Straight into corners, no , 
shear lip. 

Tension normal to rolling 
direction. 

SN-2 6 

. Yes 

Straight for short distance. 

Loading tab fracture. No 
load readout in movies. 
, 
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The fixture for holding the CSD had knife edges outside the plane 
of the specimen and bending produced some twisting of the gage that 
translated into an output. Once the crack turned so it was no 
longer at 45 deg to the loading axis, a response from the CSD was 
to be expected. 

Similar behavior was exhibited by SN-27, which was* to be a 
pure cyclic test. There should have been no response from the 

crack opening displacement (COD) gage. The COD gage was mounted 
against knife edges at the midplane of the specimen. Ideally even 
if there were bending, there should have been no response. But 
the gage holder was adhesively bonded through the entire thick- 
ness of the specimen. Defects in this bond, torsion, or bending 
in the specimen would translate into COD gage response. There 
was such response. 

If such imperfections in the testing operation were present 
as described for the five specimens discussed, it is only reason- 
able to assume that they existed in the other tests as well where 
both COD and CSD responses were expected. Some doubt exists there- 
fore as to the validity of the COD and CSD measurements for any- 
thing except as a tool for monitoring the progress of the tests 
.and qualitatively checking the response of the specimens. 

As explained in the following descriptions of test procedures, 
scribe marks were used to establish the crack growth at the end 
of each block of cyclic loads . These scribe marks were placed 
where a crack penetrated the surface. Generally it was suspected 
that the crack was longer beneath the surface because there was a 
"dimple" on the surface. When the specimens were fractured this 
suspicion was confirmed. The crack fronts were not straight 
through the specimen from surface to surface, but rather exhibited 
tunneling. Tunneling was greatest when shear lips were most 
evident and was least when there were no shear lips. 



ANALYSIS AND DATA INTERPRETATION 


Data Reduction Techniques 

Calculation of useful quantities such as stress intensity or 
crack growth rate requires considerable manipulation of primary 
data. Also, various analyses for stress intensity and related 
parameters use differing definitions for crack length, stress, 
etc. To permit verification of the interpretations that have 
been made and establish the consistency of the analyses, the 
manner in which the data was handled is reviewed below’. 


The key to calculation of R curves and strain energy release 
rates is calculation of the stress Intensity, K. In general, 

K = Co/ira [1] 

where C is a coefficient defining specimen and/or crack shape, 
and a is some measure of stress in the flaw area (how this is 
computed will also influence C) , and a is a measure of crack 
length. 


Stress, intensity at fracture, straight cracks - When the 
flaw. was not branched or curved and stayed straight during sub- 
critical crack growth under static testing, K values were com- 
puted using Holston's analysis (ref. 12). This approach was 
applied to SN-11, 15, 17, 18, 24, 25, and 26. It was also ap- 
plied to SN-2 and 6 to see if it gave different results from 
other approaches described later. For Holston's analysis. 


and 


K, 



+ P 
2A 


H 


Via 


[ 2 ] 


hi ' [3] 

2A 

where and are shown in figure 8, A is the gross area of the 

loading tab (equal to 20t for customary units or 50. 8t in SI 
units), and a is one-half the total crack length. Using a thus is 
the same as averaging the growth at each end of the crack. 


Stress intensity at fracture, branched or curved cracks - When 
the flaw was branched or curved due to previous cyclic testing, or 
grew parallel to one of the loading axes during static test, K values 
were computed using Iida's analysis (ref. 9). This approach was ap- 
plied to SN-2, 4, 5, 7, 9, 12, 14, 16, and 23. For Iida's analysis, 

K i - M SI 141 
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K II M sil A 


[5] 


where M and M are shown in figure 5, P is the load perpendic- 
bJL bli a 

ular to the predominant crack direction or branch of interest, A is 

40.23t in SI units or 15.84t in customary units, and a is one- 

half the projected crack length perpendicular to the load direction* 

The numbers 40.23 and 15.84 arise in the following manner. The Iida 

analysis for a 45° slant flaw assumes the original flaw length is 

one-half the plate width. It also assumes one ' principal stress is 

parallel to the plate length while the others are zero. To adapt 

Ihe analysis to the Maltese Cross specimen and crack configurations, 

an equivalent plate width ^ 2 b^j equal to 12 times the slant flaw 

length was used to compute Iida f s A, the ratio of projected crack 
length to plate width. For slant flaws that were nominally 5-cm 
( 2 -in.) long, the equivalent plate width was nominally 60 cm (24 
in.). For branched cracks, each branch was handled individually 
and only the load perpendicular to the branch of interest was con- 
sidered. To obtain the stress required in Iida T s computation of 
stress intensity (which for Iida is simply P/2bt) , Holston's finite- 
element analysis was used for the specimen configuration in the un- 
flawed state. The analysis yielded the following relation for 
stress at the center of the plate under uniaxial load: 

P 

Stress parallel to load direction “ a a = ^5 84t (customary units) [ 6 ] 

P 

° r Al3i (SI units) 

p 

' a 

Stress perpendicular to load direction = a, - - ~ (customary 

b 151 * 8t units) [7] 

P 

° r “ 3§576i (SI Units) 

where t is plate thickness. The accuracy of these relations was 
established with data obtained during checkout of the mixed- 
mode testing fixture. The strain measured at that time at the 
center of a 2.578-cm (1.015-in.) thick 2024 aluminum biaxial speci- 
men under a uniaxial load of 444.8 kN (100 kips) was 610 yin. /in. 
in the direction of the load. Computed strain using plane stress 
equations, a modulus of the aluminum of 7.17 x 10® N/cm 2 (10.4 x 10® 
psi) , and the above relations between load and ■ stress , is 619 x 10 6 . 
Agreement is within 1.5 %. It is interesting that the configuration 
of the specimen resulted in a small compressive stress transverse to 
the load direction under uniaxial load. For specimens SN-2 and 7, 
where cracks were curved, it was necessary to approximate the curved 
shape with straight lines as shown in figure 25 to get an equivalent 
plate width for computing X. 


47 



The 45° portions of the flaws for these two specimens, when so ap- 
proximated, were 8.15 cm (3.21 in.) for SN-2 and 9.09 cm (3.58 in.) 
for SN-7. 


Stress intensity at fracture, arc crack - The crack shape in 
SN-6 was an arc. To compute the stress intensities at fracture and 
for drawing the R curve, the analysis in reference 14 page 21.1 was 
used. Figure 28 is taken from that page. The value of stress for 
use in the equations for and K was computed from |P + P j / 2A 


where A equaled 44.91t in SI units or 17.68t in customary units. 
The numbers 44.91 and 17.68 come from Holston's finite-element an- 
alysis for the specimen configuration in the unflawed state under 
equal biaxial loads. The analysis yielded the following relation 
for stress at the center of the plate: 


Hydrostatic stress = (customary units) [8] 

P 

or 44.9it (SI units) 

where P is applied in two perpendicular directions. As above, this 
relation for stress was verified during machine checkout. Strain 
measured at that time on the 2.578-cm (1.015-in.) thick checkout 
specimen under 444.8 kN (100 kips) biaxial load was 330 to 340 x 
10“ 6. Using plane stress equations, a modulus of 7.17 x 10® kN 
(10.4 x 10® psi) , and the above relations between load and stress, 
the computed strain was 359 x 10“®. Agreement is within 8%. In 
the arc crack analysis, R is the radius of the arc. For SN— 6, the 
arc that seemed to best approximate crack shape had an R = 18.61 cm 
(7.327 in.). 


for calculating G - Note that in figures 5, 8, and 28 

the variation with a of values of the coefficient of cr/ira are dif- 
ferent for and K^. Thus, the ratio of to changes with 

a, even if the loads remain constant and the shape of the crack 
does not change. When R curves were drawn, as explained in the 
following section, according to 

R = jj ( K j 2 + K XI 2 J [9] 

for any set of corresponding loads and crack lengths, the ratio of 
experimental to was noted. When the G curve was drawn (again, 

see the following section) , this changing ratio was taken into ac- 
count in the following manner. For Holston's analysis, 

s ■ W + Ku 2 ) no, 

’ ^(h 2 + “V) mi 
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- -Hr^ K i 2 1 1 + n 1 - *i 2 "b 2 ¥a f 1 + ” 2 ) 

where n was the experimental value of K /k for each value of 
For lida's analysis, ^ 


[ 12 ] 


a. 


G = 1 E U K x 2 (l + n 2 ) « “ E - ■ M SI 2 a 2 ?ra (l + n 2 j [13] 


where again n was the experimental value of Because of the 
Xida formulation, n also equaled For the arc crack anal- 
ysis. 




ttR sin a[ 


a 

cos 2 

ll + sin 2 — 


|l + n 2 j [14] 


where, as before, n was the experimental value of For this 

formulation, n also equaled tan (a/2)* 


AG for crack growth rate curves - The formulation used for AG 
for crack growth rate curves was 


AG = G - G . 

max mm 

where G was the value of G at the peak load, and G . was the 
max r mm 

value at minimum load* 


Now 

and 

so that 


max 


mm 


1 - u" 


1 - u 2 


K. 


Imax 


+ K. 


Umax 


/K 


Imin 


+ K. 


Ilmin 


AG = 


1 - 


K. 


Imax 


+ K. 


Umax 


- K. 


Imin 


- K. 


I Imin 


• 2 ) 
m ) 


[15] 

[16] 

[17] 

[18] 


Note that this is not the same 'as 

4G ,- l^L |-( 41Ci )2 + ■ [19] 

The value used to compute K values in equation 18 was the 
average value between the beginning of a block of cycles and the end. 
Whether the actual or projected length was used depended on whether 
the Holston or Iida analysis was used, which in turn depended on the 
shape of the crack as it grew. For SN-3, 4, and 6, only the Holston 
analysis was used. For SN-5 and 8, only the Iida analysis was used. 
For SN-1, 2, 7, and 9, both were used because the crack changed 
character as it grew. Similarly, the value of stress used in the 
computation of depended on whether crack growth was being ana- 
lyzed by the Holston or Iida approach. As before, for the static 
tests, for Holston, 



[ 20 ] 


and 


o„ = 


P + P 
V H 


N 2A 


T = 


P - P 

fH 

2k 


[ 21 '] 


where A = 20t (customary units) [22] 

or A = 50. 8t (SI units) 

and, for lida, equation 6 was used. The calculation of Aa/AN was 
always the average growth of the two ends of the flaw and both 
faces of the specimen. 


Compliance - The general definition of compliance is deforma- 
tion per unit load . Compliance was computed as follows for crack 
opening displacement (COD) : 


A (COD) 

K*“h)/ 2 

and for crack shear displacement (CSD)<; 

A(CSD) 

(“v - “e )/ 2 


[23] 


[24] 


where AP was the difference between the maximum* and minimum loads 
and A(CSD) or A (COD) the difference betpeen the maximum and mini- 
mum gage response for the cycle for which compliance was being 
computed . 


Static Fracture Experiments 

A number of static fracture experiments were conducted in this 
program. In some cases, the experiment was conducted after mixed- 
mode cyclic loading with reversed shear. In others, the experiments 
were conducted immediately after the original fatigue sharpening of 
the crack. As described earlier, an attempt was made in all experi- 
ments to record crack length and load on both axes as functions; of 
time throughout the final phase of the • experiment . Because the 
thickness of materials ' tested was such that plane strain conditions 
for fracture were not maintained, significant subcritical crack 
growth was evident before final failure occurred. This implies the 
presence of a significant amount of plastic deformation in the frac- 
ture process, so it was decided to analyze the failure of these 
specimens in terms of the R-curve concept. 

Figure 29 is a plot of G, strain energy release ’rate, versus a, 
crack length. The line labeled R Curve represents the amount of 
energy per unit crack area required to drive the crack through the 
material. This curve conceptually includes the energy necessary to 
form the large plastic zone associated with fracture in the speci- 
mens tested in this program. When high-speed motion— picture data 
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are available on crack length and load as the crack grows toward 
failure, it is possible to construct the R curve for each experi- 
ment as follows: 

1) The point labeled a^ is initial crack length; 

2) The line labeled o ^ represents a theoretical calculation of 

G as a function of crack length for stress level The 

stress level and corresponding crack length are known 

from high-speed motion-picture data, so it is possible to 
compute a value for G that locates the data point on the 
R curve; 

3) This process is repeated for stress levels etc > un ~ 

til a point of tangency is found between the. R curve and 
the equation for G as a function of a; 

4) The tangency point is taken to be the critical point at 
which rapid crack propogation occurs* This point identi- 
fies the critical stress, a , and critical crack length, 

a ; 
cr 

5) These quantities are used to compute G , the critical 
strain energy release rate. 

To construct the R curve from experimental data as ^described 
above, a method to calculate G is necessary. In the previous sec- 
tion, methods for calculating stress intensity factors K^. and 

and the strain energy release rate were described. To repeat, in 
constructing the R curves for each static fracture experiment, the 
following equation was used: 

G - (V + V) '“l 

This expression correctly represents G as long as K^. and values 

are computed for crack growth in the direction in which the crack 
actually is growing. As described in the previous section, every 
attempt was made to do this. 

Figures 30 through 45 present R curve constructions for speci- 
mens on which there were enough data. Figures 30 and 31 present R 
curves for sample SN-2 calculated by two different analysis proce- 
dures. Figure 30 used the curves of Iida, and figure 31 used the 
calculations of Holston for the same specimen. Little difference 
is noted between the two R curves produced. Also, the critical 
value of G and critical crack length predicted by the two analysis 
procedures are very nearly equal. This lends credence to the use 
of Iida’s analysis, which is valid only for uniaxial loading in a 
situation in which loading is actually biaxial. The results in 
figure 36 were computed using and expressions for a two- 
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Crack extension force, G, and 
crack extension resistance, R 





'Figure 30* - R curve for SN-2, lida analysis. 









Crack extension force, G, and 
crack extension resistance, R 



Figure 32. - R curve for SN-4. 
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Figure 34. - R curve for SN-18, 



Crack extension force, G, and 
crack extension resistance, R 



Figure 35. ~ R curve for SN-5. 




Crack extension force, G, and 
crack extension resistance, R 



Figure 36, - R curve for SN-6, arc-crack analysis. 






dimensional arc crack. Figure 37 analyzes the same experiment using 
Holston T s analysis. Again, differences between predicted critical 
G and critical crack length are small. Figures 39 and 40 present 
two analyses’ of sample SN-12. In this case, there is considerable 
difference between critical G and predicted critical crack length, 
but it is important to note that this sample was loaded in pure 
shear, meaning that the vertical load was compressive. Sample SN-12 
failed by buckling rather than by crack propagation, making the R 
curves in figures 39 and 40 subject to question. Each of the re- 
maining R curves was computed based on a single analysis described 
in the previous section. 

R curve results are summarized in table VI, where they are 
categorized according to the alloy, macerial thickness, direction 
of crack growth relative to final rolling direction, and the ratio 
of mode I to mode II stress intensity that actually occurred in 
the experiment. In some experiments, data necessary to construct 
an R curve were not available. In these cases, values of critical 
G listed in table VI were computed from an estimate of the load at 
fracture. 

Specimens SN-7 and 9 failed in rapid fracture while cyclic 
stress was being applied. In these cases, the load at fracture 
was recorded on the last cycle before failure, and the crack length 
at failure was estimated from the last recorded crack lengths be- 
fore failure. Specimens SN-23 and 26 were rising-load tests and 
showed little sub critical crack growth. In both cases, the load 
at fracture was recorded. For SN-23, three frames of motion pic- 
ture data are available from which critical crack length was esti- 
mated. The loading tab broke off SN-26 before fracture, so the 
estimated G value is well below the appropriate critical G. 

Table VI gives values for net section normal and shear stress. 
Several experiments were conducted with the fracture direction 45° 
from the final rolling direction, and the loading was pure mode II, 
or the ratio was approximately 0.4. In these tests, the 

failure was nearly flat and in the original plane of the crack. 

Note that net section shear stress is well above the expected shear- 
ing yield stress of approximately 0.6 times the tensile yield stress, 
so it must be concluded that these specimens (SN-11, 17, 18, 24, 25, 
26) failed by net section yielding. Note that in some cases there 
is wide variation in the value of G ^ within a given alloy/ thickness 

category. For a given alloy, it is expected that the following par- 
ameters should affect the critical value of G: 

1) Material thickness - Material thickness influences the con- 
straint at the crack tip and thus the amount of plastic de- 
formation that can take place during the fracture process; 

2) Crack growth direction - The direction in which the crack 
propagates relative to the rolling direction of the material 
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Figure 40. - R curve for SN-12, Holstcn analysis 




TABLE VI 


R CURVE SUMMARY 


Specimen Material 
no. 


Crack growth 
direction re- 
lated to 
j final rolling 
direction 


kH-cm 1 ^ 2 / | ksi-m .*/ 2 | kN I ksi-in*^ 


2219 

Parallel 

t=0. 64cm 

Either 

or 

45 deg 

t^o.25 in. 

Perpendicular 


Parallel 
t=1.55 cm | Parallel 
or 

t-0. 61 in 


Parallel 

57.14 

Parallel 

51.10 

Perpendicular 

102.38 

45 deg 

72.85 

Perpendicular * 

76.92 

Parallel 

47.47 


22.42 

20.00 

95.61 

95.21 

>68.14 


20.40 

18.20 

87.01 

86.65 

>62.01 


Specimen 




G 

cr 

Nomin 

kN-cm 

kip-in. 
in. ^ 

Net section 
Normal stress 



kN/ cm 2 

ksi 

0.88 

0.50 


0.96 

0.55 

NA 

NA 

0.86 

•0.49 

NA 

16.3 

0.46 

0.26 

11.2 

NA 

0.35 

0.20 

NA 

NA 

0.37 

0.21 

NA 

NA 

0.'68 

0.39 

10.6 

15.4 

0.89 

0.51 

=0 

-0 

0.40 

0.23 

NA 

NA 

0.33 

0.19 

NA 

NA 

1.33 

0.76 

NA 

NA 

0.75 

0.43 

2.81 

40.7 

0.81 

0.46 

NA 

NA 

0.33 

0.19 

NA 

NA 

1.16 

0.66 

=0 

-0 

1.14 

0.65 

-0 

-0 

>0.68 

>0.39 

11.2 

16.2 


Net section 


NA 

NA 


Yield , 
Stress 


NA 

NA 


NA 

NA 

NA 

=0 

NA 

NA 


ksx 


ksi 

NA 

38 

55 

NA 

38 

55 

39.6 

38 

55 

NA 

38 

55 

NA 

38 

55 

NA 

38 

55 

37.1 

38 

55 

42.6 

38 

55 

NA 

45 

65 

NA 

45 

65 

NA 

45 

65 

=0 

45 

65 

NA 

45 

65 

NA 

45 

65 

48.3 

45 

65 

47.8 

45 

65 

36.9 

45 

65 


55 | Arc crack 


55 Buckled 


_ 37.7 
28.5 
54.8 
Pure K-j. 

3.4 

2.4 

Pure K t 


65 No R curve 
65 No R curve 


65 No R curve 


65 No R curve 


6.65 2.62 
6.05 2.38 
6.48 2.55 
7.87 3.10 


13.11* 5.16* 

11.07* 4.36* 

4.67 1.84 

6.05 2.38 

5.08 2.00 
4.42 


*Extrapolated. 


[*!•] 









































Crack extension force, G, and 
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Figure 45. - R curve for SN-25* 





might be expected to influence the critical value of G be- 
cause of differences in material properties with respect to 
direction relative to rolling direction; 

3) Ratio of Kj. to - This ra'tio should have some effect on 

the critical value of G because it is expected that develop- 
ment of plasticity is different in mode I than in mode II. 


To gain some physical insight into the effects of plasticity in 
mode I and mode II, it is instructive to examine the plastic zone 
dimensions predicted by use of Dugdale-type strip models. The fol- 
lowing expressions taken from reference 12 represent the estimated 
plastic zone dimensions for a condition of pure modes I and II: 


h ■ 1 (W) 2 ' 1251 

L II = 2 ( K Il/°ys) 2 ^ 26 -* 


where L refers to the plastic zone dimension. Note that, for a given 
level of K for pure mode I or II, the plastic* zone dimension for pure 
mode II is estimated to be four times greater than that for pure mode 

I. This implies that, there should be greater plastic energy to over- 
come in initiating fracture for pure mode II than for pure mode I. 

One therefore concludes that, as the loading condition traverses from 
pure mode I through conditions of mixed-mode loading toward pure mode 

II, the energy of fracture, or G , might be expected to increase. 

cr 


IThe effect of thickness on critical values of G has been long 
recognized for conditions of pure mode I. As material thickness de- 
creases, f’he amount of constraint at the crack tip also decreases, 
thus allowing plastic deformation to increase. The result is an in- 
crease in the critical value of G as thickness decreases. 

Figures 46, 47, and 48 are plots of the critical value of G de- 
termined from the R curve plotted against the angle 8 defined by 

0 = arctan (K^/K^) [271 

where the subscript u indicates fracture. 

In these figures, data points are annotated to indicate the di- 
rection of crack propagation relative to the final rolling direction 
of the specimen being tested. Figure 46 presents the result for 
tests of 1.55-cm (0.61- in.) thick 2219. Note ‘that the variation 
in the critical value of \G with Kj/Kjj is in the direction predicted 

above, and there appears to be little separation of the data with 

respect to the crack growth directions present. In figure 47, the 

results are presented for 0.64-cm (0.25-in.) thick 2219. In this 

case, the variation seen in figure 46 is not evident. The G shorn 

cr 

was computed assuming flat cracks. However, examination of the* 
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Figure 46. - Effect of mixed-mode ratio 
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Figure 47. - 
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Figure 48. - Effect of mixed-mode ratio, 7075 1.27-cra (0.50-in.) thick. 




as determined from the R curves .in 


specimens in the area of a 
r cr 

figures 35 and 36 showed that the cracks had twisted approximately 
45°. A better estimate of G cr for SN-5 and 6 taking account of the 

rotated crack and resulting presence of mode III. would be about 60% 
of the values shown in figure 47. Since the interaction diagrams 
in the next section involve only and , consideration of 

and crack rotation were ignored. Furthermore, specimens SN-5 and 6 
were the only ones exhibiting this behavior. Figure 48 is a simi- 
lar plot for the 1.27-cm (0.50- in.) thick 7075. In this case, the 
proper trend with respect to angle 6 is seen for tests in which 
crack growth direction was 45° to the final rolling direction. 
However, for cases in which crack growth was perpendicular to the 
final rolling direction, the variation is not consistent with the 
model described earlier. Also, data for the parallel propagation 
direction are nearly horizontal and do not support or contradict 
the proposed model. In figure 48, note the legend corresponding to 
the data point for specimen SN-26. This specimen did not undergo a 
rapid fracture; rather, the tab broke off the specimen at a G level 
of about 0.7 kN-cm/cm 2 as indicated in the figure. The data point 
is included in a way that indicates that G for that experiment would 
have been a value greater than 0.7. 

One can conclude from these results that the critical value 
of G for pure mode II loading situations generally exceeds the 
critical value of G for pure mode I loading situations; and that 
the critical value of G for intermediate conditions of mixed-mode 
I and II loading lies between the two extremes. This statement 
is generally supported by the data; however, the 7075 data for 
perpendicular crack growth direction are not in agreement and re- 
sults shown in figures 47 and 48 for parallel crack growth di- 
rection are inconclusive with respect to this point. 

The R curves described in this section could have been plot- 
ted in terms of S in the strain energy density theory instead of 
G as used here. This was done for one specimen, SN— 17, and the 
plots are shown in figure 33. Since there was almost no difference 

in the determination of a , and because there was little differ- 

cr 

ence to be expected theoretically, as described in the state-of- 
the-art review, analyses using S were not done. 


Interaction Diagrams 

Figures 49, 50, and 51 are plots of all static fracture data 
obtained in this prpgram in the format of an interaction diagram. 
There is significant scatter in the data when plotted this way. 
However, if one separates these data according to crack growth di- 
rection relative to the final rolling direction of the specimen, 
one finds a tendency for the data to be stratified with respect 









to crack growth direction. To normalize the data, several critical • 
stress intensity values from these figures were selected to repre- 
sent critical K values for each group of material thicknesses ..and 
crack growth directions relative to final rolling direction. Critic 
cal stress intensity values selected are summarized in Table VII. 
Table VIII presents a series of data from reference 15 for the 

same materials in different thicknesses. Note that numbers se- 
lected here are consistent with other data on these materials. 

Figure 52 is a normalized plot of data that were first presented 
in figures 49, 50, and 51. Data in these figures were categorized 
according to crack growth direction and were normalized using the 
appropriate K^ c value taken from Table VII, In figure 52, note that 

normalizing in this way significantly reduces data scatter. Also 
note that there is a significant degree of interaction evident in 
these data, and the interaction is nearly as severe as the straight- 
line interaction noted by Shah (ref. 1) and Liu (ref. 4). Further 
note that the critical value if is significantly greater than 

that of Kj. Figure 52 includes plots of predicted failure envelopes 

according to the various failure theories described in the state-of- 
the-art section. Values from the data show greater capability for 
combined mode stress intensities than the theories predicted for high 
proportions of shear stress intensity. However, in this region, 
net section stress exceeded yield stress in shear, as noted before. 
For small proportions of shear stress intensity, there is good cor- 
relation between data and theories. 

Figure 53 is an interaction diagram in which the axis has 

been normalized according to the appropriate values of taken 

from Table VII, and the K^-axis has been normalized according to 

the values of taken from Table VII. A further dcrease in 

data scatter is seen but it is somewhat artificial because of the 
high degree of normalization. 

The critical shear intensity, tends to be greather than 

Kj c for the series of experiments conducted, and there is signif- 
icant interaction between modes I and II. Further, it was found 
that the various theories of fracture correlate well with the 
data for low proportions of shear but underestimate stress in- 
tensities at fracture for high proportions of shear. It must be 
recognized that, in these experiments, a significant degree of 
plasticity influences the fracture process under high proportions 
of shear, and that the theories of fracture described only apply 
to conditions of brittle fracture. 
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TABLE VII. - PURE MODE CRITICAL STRESS INTENSITIES, THIS PROGRAM 


Material 

Thickness 

K t 

Ic 

K tt 

lie 


cm 

in. 

Parallel 

45 deg 

Perpendicular 

kN-cm 1 / 2 

ksi-in. */ 2 




kN-cm 1 / 2 

cm 7 

ksi- 
in. r / 2 

kN-cm 1 / 2 

cm z 

ksi- 
in . 1 / 2 

kN-cm 1 / 2 

cm z 

ksi- 

in. 1 / 2 

cm 2 


2219- T87 

0.64 

0.25 

85.7 

78 

85.7 

78 

NA 

NA 

119.8* 

109* 

2219-T87 

1.55 

0.61 

53.8 

49 

53.8 

49 

NA 

NA 

85.7 

78 

7075-T7351 

1.27 

0.50 

51.6 

47 

72.5 

66 

89.0 

81 

95.6 

87 


^Estimated 

NA = not available. 


TABLE VIII. - VALUES FROM DAMAGE-TOLERANT DESIGN 
HANDBOOK (ref. 15) 


Material 

— 

Thickness 

K 

c 

cm 

in. 

kN-cm 1 / 2 

cm 2 

ksi-in. 1 / 2 

2219-T87 

0.152 

0.060 

84.6 - 86.8 

77-79 

2219-T87 

0.254 

0.100 

82.4 - 115.4 

75-105 

2219-T87 

3.175 

1.250 

39.6 - 49.4 

36-45 

7075- T6 

0.102 

0.040 

57.1 - 75.8 

52-69 

7075-T6 

0.127 

0.050 

63.7 - 78.0 

58-71 

7075-T651 

3.505 

1.380 

22.0 - 33.0 

20-30 













Fatigue Experiments 


A number of fatigure studies were conducted for which crack 
length was recorded as a function of number of cycles and loading. 
Figure 24 presented diagrams that describe variation with time of 
vertical and horizontal loads applied to cracked specimens. The 
first part of each cycle in figure 24 is referred to as part0; 
the second is referred to as part (2) . In these fatigue experiments, 
the initial precrack was oriented 45° to the vertical and horizontal 
axes. However, in most tests, crack branching occurred as soon as 
mixed-mode cycling with reversed shear commenced. Branching occur- 
red in such a way that one branch grew predominantly normal to the 
vertical loading; the other grew predominantly normal to the hori- 
zontal loading. This meant that, on a given crack tip, if ap- 
plied loading produced predominantely mode I with some mode II 
on one part of the cycle, then almost no stress intensity of 
either kind was produced during the other part of the cycle. 


Because there was some degree of mixed-mode loading, fatigue 
data were analyzed in the following way. It was decided to charac- 
terize crack growth in terms of AG instead of the usual AK. AG quan- 
tities were computed for parts (T) and (2) of the cycle using the fol- 
lowing expressions: 

AG ® ,(D = [( K I 2 + K II 2 )max " ( K I 2 + K II 2 )min]0 ,0 [28] 

' i 

Next, it was assumed that a fatigue law of the type 

da/dN = C f ( AK x ) n • [29] 

could be extended to mixed-mode situations by replacing the term AK T 


with the equivalent relationship in terms of AG given by 

\n/2 


(ag) ” /2 


[30] 


Next, it was considered that crack growth occurs during both parts of 
the cycle plotted in figure 24. Accordingly, crack growth rate is 
assumed to contain contribution due to both parts (T) and (2) of the 
cycle and is expressed as 

da/dN = C (~-^) J AG Q) n/2 + ( AG 0) n/2 ] 


[31] 


which may be rewritten- in the form 

i n/2 


da/dN = C 


(r HT 


AG 


.n/2 


(D 


1 + 


/AG 

AG 


_ vn/2- 

0/ J 


[32] 



The next step in the analysis is to recognize that the ratio of AG 
to is a small number and may be neglected compared to unity. 


© 


Based on the foregoing analyses and series of assumptions, it 
is concluded that the mixed-mode fatigue data obtained in this pro- 
gram should correlate with an expression having the form 


da/dN = C 



[33] 


where parameters n and C may be taken from mode I fatigue data ob- 
tained in previous experiments. The foregoing discussion presents a 
way in which crack growth data may be interpreted in terms of AG for 
mixed-mode situations. Equations of a type other than that of equa- 
tion 29 could have been used as a basis for the discussion as well, 
with replacement for AK made in terms of AG in a similar manner. 

For example, in the following discussion, some comparisons are made 
with Forman’s (ref. 16) equation presented in terms of AG instead of 
AK. 


Figures 54 through 59 present fatigue data from this study for 
2219 aluminum in the form of plots of AG versus Aa/AN. The figures 
also show plots of Forman’s equation evaluated for 2219 aluminum in 
mode I fatigue tests. Data from this study involving fatigue crack 
growth under mixed-mode loading conditions were found to agree very 
well with the Forman equation when it is interpreted in terms of AG 
as described above. 

Figure 60 through 63 are plots of AG versus Aa/AN for 7075 alum- 
inum. These data are compared to Colliprist data (ref. 17)* for the 
same material subjected to mode I fatigue loading. Very good agee- 
ment is found between the mixed-mode fatigue data from this program 
so analyzed and data based on Collipriest ’ s mode I fatigue stud- 
ies for the slower growth rates. 

In the upper ranges of Aa/AN, data from this program con- 
sistently fell below Collipriest ' s pure mode I data. The dif- 
ference cannot be attributed to the presence of mode II effects 
because there was very little mode II loading in the latter stages 
of the tests. It is more likely attributable to differences in 
specimen dimensions and shape and/or stress levels relative to 
yielding under which crack growth was fostered. 

It is also interesting that in figure 63 the growth rate for 
SN-27 was always faster than predicted by Collipriest. In SN-27, 
the crack was grown under essentially pure conditions except 

*Collipriest’s curve was not used. A curve was fitted through 
Collipriest 1 s data points in the range Aa/AN covered by the 
tests on this program. 
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that a compressive load applied parallel to the crack growth direc- 
tion on every cycle was equal in magnitude to the tensile load. 
Apparently this compression had considerable effect on growth rate, 
although no theoretical treatment takes such loading into account. 

It is concluded from these results that the method presented 
here for analyzing mixed-mode fatigue data is reasonably correct for 
the experiments conducted. However, it is also recognized that the 
amount of mode II loading present in these fatigure experiments 
should not be expected to have a highly significant effect on the 
results . 
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Figure 56. - Crack 
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Figure 57* - Crack growth 
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Compliance 


Figures 64 through 73 show compliance data. Both analytical and 
experimental results are given for comparison. These results were 
not used in any other calculation; crack lengths were measured di- 
rectly, and stress intensities were calculated using measured loads 
and crack lengths. Compliance data were used to monitor tests and 
provided a direct comparison between test and analysis. 

Test data points shown were obtained by dividing measured COD 
and CSD by + P R j and - P^j respectively, except in cases 

of single-mode loading. In. cases of pure mode I loading (SN-1, 2, 

6, 7), P = P , and CSD should be zero as long as the specimen is 
V H 

symmetric. Thus, shearing compliance, 2(CSD)y^P v - P^j , is inde- 
terminate (0/0). For these cases, the ratio CSD/COD is shown along 
with opening compliance. Deviation of this ratio from zero indi- 
cates asymmetry of the specimen and/or deviation of loading from 

pure P T , = P (see section on experiment observations and primary 
V Jti 

data). For pure mode II (SN-27), P - -P , and COD should be zero. 

V H 

The ratio COD/CSD is shown along with shearing compliance. Ana- 
lytical results were obtained from the following expressions (see 
Appendix A for development) : 

COD - 4 ^ ~ u ' a + c H )v [34] 

CSD - 4 - k ; * H(o v - ° H )u [35] 

with = P^y/20t and = P^20t (customary units) 

or a v = P^SO.St and - P^y^jO.St (SI units) 

where t denotes thickness and 50.8 or 20 is the length used in nor- 
malizing analytical results. Theoretical compliances are shown by 
heavy lines, and the extent of the lines indicates range of appli- 
cability. Analytical results are shorn where the crack remained 
straight or branched with equal branches. For other configurations 
(2 , unequal branches*, or curved cracks), analytical expressions for 
stress intensities are available, but not for crack displacements. 

These figures show general agreement between tests and analyses. 
The correlation deteriorated as the cracks grew into configurations 
different from those idealized in the analyses. In cases of single- 
mode loading, pure K^., or pure the other crack displacement 

should have been zero. The figures show some small values. As ex- 
plained above, these may have been due to asymmetry or bending of 
the specimen, small differences in the loads on each axis, and/or 
possible gage-holder mounting problems. 
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Figure 73. - Compliance for SN-27 





CONCLUSIONS 


To develop crack-growth and fracture toughness data under mix- 
ed mode conditions, an experimental investigation was undertaken 
involving, in the main, room- temperature cyclic and static tests of 
flawed 91x91-cm (36x36- in.) aluminum plates. General conclusions 
that we believe are supported by the test results are set forth in 
this section. 

During cyclic testing, it was impossible to maintain a high 
proportion of shear-mode deformation on the cracks. The cracks al- 
ways branched when the shear was reversed on each cycle, resulting 
in crack fronts with very little When the shear was not re- 

versed on each cycle, the cracks simply grew in an orientation that 
again resulted in almost no at the tip. When there was no shear, 

the cracks turned to follow the final rolling direction of the mater- 
ial, When cracks branched, the branch parallel to the final rolling 
direction eventually grew faster. 

Mixed-mode crack- growth* rate data were analyzed in terms of 
energy release rate, G, taken as a function of both K^and , ver- 
sus growth rate. Test results compare reasonably well with single- 
mode data from Forman and Collipriest, when their data are handled 
in the same fashion, keeping in mind that the shear component of G 
was relatively small. There was no noticeable difference in the 
effectiveness of the analyses used on crack growth data. Even though 
straight and branched cracks were analyzed using Hols ton factors, 
while Z-shaped cracks were analyzed using lida factors, there was 
not an abrupt change in computed quantities derived from the test 
data when the cracks changed configuration and the analytical ap- 
proach was correspondingly changed. In fact, in transition ranges, 
for example when cracks turned, the various approaches gave nearly 
the same computed values of K or G. 

Measured crack opening and crack shear displacement agreed well 
with the straight and branched crack analyses over the range of ap- 
plicability, that is, while the cracks were straight or the branches 
grew equally, and when significant displacements were expected from 
the load scheme applied; However, there was measurable opening and 
shear displacements when theoretically they should have been zero, 
suggesting that there were deviations from ideal conditions of speci- 
men shape, load application, and gage mounting. 

Under static loading to fracture, straight 45° cracks remained 
straight when the shear stress intensity exceeded the normal stress 
intensity at the crack tip. Thus, we were able to obtain fracture 
under pure shear. When applied shear was less than the normal 



stress, the crack grew, or turned to become, perpendicular to the 
direction of higher load. When the loads on each axis were equal, 
the crack turned parallel to the final rolling direction. If the 
crack was branched before the static test, the branch perpendicular 
to the higher load was the one that grew to instability. When sta- 
tic loading of equal magnitude on each axis was applied to a curved 
crack, the crack grew to instability without further change in di- 
rection. 

It was shown theoretically that there is little difference in 
the prediction of fracture behavior between the maximum tangential 
stress theory minimum-strain energy density theory (S) , and 

an extension to the theory based on energy release rate (G) . There- 
fore, all data were reduced in terms of G because it was simpler. 

In the range near pure K^, the experimental G values fell both in- 
side and outside the theories. Values of G for pure K_ were a\>- 

cr r II r 

proximately 50% higher than G £r for pure K^, which is contrary to the 

three theories mentioned above. However, the net section stresses 
were well above the shear yield stress when fracture occurred with 
high shear. For the 7075 alloy, values of G were also affected by 

the relation between flaw growth direction and final rolling direc- 
tion of the material. This dependence was not found in the 2219 al- 
loy. For both, there was very strong interaction between the normal 
and shear mode stress intensities, with a large reduction in K^. re- 
sulting from the application of K__. Again however, the conclusion 

must be tempered by the fact that fracture under conditions of high 
occurred when the net section stresses were well within the in- 
elastic region. 

Choice of the Maltese cross specimen shape for this investiga- 
tion served the program well. This configuration permits any com- 
bination of stress intensities from pure through pure with 

a 45° flaw orientation. The straight-crack analysis of reference 
10 showed that stress intensities are nearly constant for a range 
of crack lengths. Branched crack results given in Appendix A show 
similar behavior. These analytical results were obtained via a 
special crack-tip finite element that properly accounts for both 
singularities. Resulting stress intensities are applicable for 
any thickness. Young's modulus, and loading combination for geo- 
metrically similar specimens. 

In addition, the checkout procedure on the test machine re- 
vealed that the area of strain uniformity (±10%) under biaxial load 
in the center of an unflawed specimen extended for about five times 
the original 5-cm (2- in.) crack length, so that flaws grew into re- 
gions of constant strain. 
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APPENDIX A 


ANALYSIS OF TEST SPECIMENS 
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I. INTRODUCTION 


The test specimens shown in Figure 1 x*ere analyzed to determine stress intensity 
factors and crack displacements. Stress intensity factors are used in studying 
flaw growth and fracture. * Crack displacements may be compared with experimental 
measured values to provide a correlation between the analysis and test. The 
analysis* considered specimen plane geometry and thickness, material elastic properties, 
flaw length, and combined loads. The main results are presented in the form of 
graphs for normalized stress intensity factors and crack displacements versus flaw 
length.. 

In performing the analysis, the finite element method was used to account for 
the irregular shape of the specimen boundary. A conventional element was used for 
most of the model and a special element used to model the region around the crack 
tip. This special element is required to properly treat the singularities occurring 
at the crack tip. The effects of flaw length were treated by making different 
finite element models for different flaw lengths. 

Several features of the analysis method and specimen geometry permitted 
simplification of the modeling and generalization of the results. Symmetry of the 
specimen permitted modeling only one fourth of it. Linearity of the analysis with 
respect to specimen thickness, Young 1 s modulus,, and loads permitted scaling of 
these parameters. Finally, through the use of superposition any combination of 
loads may be treated from two fundamental cases. These features are discussed 
further in subsequent sections. 


PRECEDING PAGE' BLANK NOT HLM^ 
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The main results from this analysis are given in the form of graphs showing 
normalized stress intensity factors and crack displacements versus flaw length* 

These results are applicable to other geometrically similar specimens for any 
thickness. Young 1 s modulus, and all load combination* 

II* SCALING 

Before proceeding with scaling, consider the relative stiffnesses of the 
grip plates and test specimens. The grip plates will be used repeatedly and 
they must remain in tolerance with an infinite fatigue life. Thus they will be 
designed to low stress levels and will be very stiff relative to the test specimen 
Consequently, the grips were not modeled in the finite element analysis and 
a uhiform displacement x*as applied to the specimen to simulate the grip plates. 

The load associated with the uniform displacement was determined from the analysis 

Two types of elements were used in the finite element models. Most of the 

» 

specimen was modeled with constant strain triangular elements. This element is 
given in many text books on finite element analysis, such as Reference 1, and its 
stiffness matrix is linear in element thickness and Young's modulus, 

| K' J = h E I f (x^ y i> V >V (1) 

where h, E-, ']/ are thickness, Young's modulus, and Poisson’s ratio, respectively, 

x and y are node point coordinates, and 
i 1 
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1/ - y for plane strain, 


( 2 ) 


V - 


1/ 

1 +v 


for plane stress. 


The region around the crack tip was modeled with a special finite element. 
Reference 2. This element is based on Williams eigenfunction expansion of an 
airy type stress function about the crack tip. It is circular in shape, applica- 
ble to Kj. (opening), (in plane sliding), and .modes of deformation, 

and contains the proper singularities. Generalized coordinates were used in the 
development and higher order terms from the expansions included. Compatibility of 

displacements at node points of adjacent conventional elements was ensured -by 

* 

using the Lagrange multiplier method in a constrained minimization. The stiffness 
matrix for this element is linear in Young’s modulus- and thickness also 

[Kj = h E [f{7,V )] (3) 

where the element geometrical parameter (element radius, r^) is contained in the 
generalized coordinates. Stress intensity factors are given by 


K i - - 2 


3/2 


E 


2(1 +y 


) /r 


*11 


= 2 


3/2 


E 


2.(1 +y 




(4) 


where and are generalized coordinates. Now consider the question of scaling 
of stress intensity factors with Young's modulus and thickness. Suppose two finite 
element models were assembled from the above elements. The equations for these models 
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would be 



( 5 ) 


where the subscripts denote the two models* Now each model stiffness matrix is 
linear in Young’s modulus and thickness, since it was assembled from linear 
elements. Thus if the models have the same plane geometry and restraint, then their 
stiffness matrices are related as follows 



If the node point loads are proportional, then 



The solution for the second model is 



and (6) shows that 

Substituting (9) and (7) into (8) gives 



( 6 ) 


(7) 


( 8 ) 


(9) 


( 10 ) 


Equations 4 show that stress intensity factors for each model are proportional to 
Young's modulus and generalized displacements. • Combining this with (10) gives 



Thus, stress intensity factors for all test specimens may be obtained by scaling 
results from a single model. In scaling, Young's modulus, specimen thickness and loads 
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are scaled. All other features must be equal. These Include plane geometry 
(specimen including flaw length), Poisson's ratio, and type of plane analysis 
(plane stress or plane strain) . 

Scaling relationships for stresses and strains throughout the specimen 
may be developed. Strains are proportional to node point displacements 
and stresses are linear in Young's modulus and strains. Combining these with 
(10) gives 



( 12 ) 


III. SYMMETRY AND SUPERPOSITION 


Both the straight crack and branched crack configurations (Figure 1) have two 
axes of symmetry^ along the diagonals of the specimens. Thus only one-fourth of 
the specimens need be modeled. Results for the remaining three- fourths of the 
specimens are obtained from symmetry relationships and results computed in the 
quarter model. These symmetry relationships are developed in this section. 

Figure 2 shows the coordinate system. P(r,9) is a generic point and 

1 

P (r,04-TT)is its image point. The structure is periodic in the angular coordinate 
9 with period It • If the loading is also periodic, then the response will be 
periodic. 


JJ r (r,0+[ t) =// r (r,0) 

.U 0 (r,9+’iT) =£/ 0 (r,e) 


(13) 


This reduces the modeling required by %. 
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Figure A-l. 



Test specimens. 



/ 


Figure A-2. — Polar coordinate system. 
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The structure is also symmetric about the line 0 = Ojjf Thus for symmetric 

loads, the response will be symmetric and for antisymmetric loads, it will be 

antisymmetric* Now any periodic loading can be decomposed into its symmetric 
( 

and antisymmetric parts as shown in Figure 3* These two cases, symmetric and 
antisymmetric, may be treated separately and the results superposed since the 
analysis is linear in loadings. This provides another % reduction in the modeling. 
However, two different sets of boundary conditions on the quarter model must be 
treated corresponding to the symmetric and antisymmetric modes of deformation. 
Development of the required symmetry/ antisymmetry relationships and boundary 
conditions is facilitated, by considering several special load cases as follows. 

For the first case consider the loading as shown in Figure 4 and denote this 
as case M a IT . Let Q(r,-0) denote the image of P(r,0) with respect to the line 
9 = 0 and assume the displacements to be positive as shown* Now consider a load 
case obtained by folding case fl a ,T along the line 9 = 0,’TT^ as shown in Figure 5^ and 
denote it case n b n . Since the structure is symmetric with respect to this line, 
the displacements for case rr b n are those of case ,f a IT folded about the line as 
shown. Thus 

u T v\-u t i Q) a -" e W b = -V Q) a 

(14) 

" r (Q >b =U r (p >a 

The third case is the superposition of cases "a" and fT b n as shown in Figure 6. 

\ 

Displacements for this case- are obtained by superimposing those of cases n a" and ,f b M , 
thus, 
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P(-Q+Tr)= P(-6) 


U r =U r (Q) a 

* u e* / 

P A- -Ug(Q) a < 


Figure A-5. - Load case b. 
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p(-e)=f(e> 




(15) 


Xj r ( p ) a+ b - U r (P >a +U r (P) b 

A. lr (Q) i+b - U r (Q) a +U r (Q) b 
e^a+b = -1' 9 + l J lj l ' P \ 
Xt 9 ( ®a+b = *•' e W) a + *V®b 


Substituting (14) into (15) gives 

14 r (P) a+b =Z( r (P) a +i/ r «) a 

^ r <Q) a+ b-X( r «)a + V P) a 

H 9< P:> a+b = ''V P) « --<V Q) a 

0 ^ a+b' = ^ 0 ^ a " M 0 ^ a 
Equations 15' show that 

r^a+b = ^ r^a+b 

XA 0^a+b = -t/ 6^^a+b 


(16) 


(17) 


as is .shown in Figure 6. Substitution of (13) into (17) gives 

■«y) rt ■( ' r (pl >a+b 

(18) 

■t- { e (Q ^a+b = “^e (P ^a+b 

]_ f 

as shown. Now P(r,9) is a generic point and the other points, P ,Q and Q, are images. 
Thus letting Q=0 in (17) and (18) defines conditions along the line 9=0,77'. Similarily, 
setting 0 ='7T/2 defines conditions along 0 =7/ /2, -'T/ll as shown in Figure 6. . For 
the straight crack configuration, this loading gives pure opening stress intensity 
(K^.) . In the branched crack configuration, both modes and are present, with 
their proportion dependent on branching angle and length. 
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The final load case is the superposition of case tT a u with the negative of 
case !t b n as shown in Figure 7. Thus superimposing displacements of case TI a IT with 
the negative of case n b tr gives 

>U r «) a . b -JU t «» a vU r (Q) b 
sj. e < p ) a _ b m -U e < p ) a e < p ) b 
kx Q (® a _ b = U e (Q) a -U 0 (® b 

Substituting (14) into (19) gives 

JU r «> a -b=U r (© a -U r ( p ) a 

( 20 ) 

-u e ( p > a - b -ll'e (p) a +l -'e<®a 
U e (Q) e-b = ^>e (Q) a + U 6 < P >a 

Comparing equations 20 shows that 

■U r «) a . b - -J.l r ®>a-b 

( 21 ) 

d- 9 a-b “ ^ e ^ a-b 

as shown in Figure 7. Setting 6 = 0 andf^/2 give conditions along the diagonals 
as shown. For the straight crack configuration, this loading gives pure shearing 
stress intensity (K ) . The branched crack gives a mixed mode (Kj. and K^) with 
different proportion than the previous case, in general. 


126 




Figure A-7. - 



i case a - b. 
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Thus a % model of the specimen along with superposition and scaling will 
provide answers for the complete specimen, all thicknesses and Young 1 s moduli, 
under any load combination. 

Now consider the question of rigid body motion and restraint. In the pro- 
ceeding discussion, it was implicitly assumed that no rigid body motion was 
present. Three rigid body modes are possible in a plane problem; two translations 
and a rotation. For the case "a + b M , sufficient restraint is provided by the 
symmetry conditions to prevent all three rigid body motions. as is shown in 

t 

Figure 6. In the case n a - b" symmetry conditions prevent both translations but a 
rigid body rotation about the origin would be possibleas shown in Figure 7. This 
rotation may be prevented by restraining an additional node point or within the 
crack tip element as it contains all three rigid body modes explicitly. In this 
study a uniform displacement was imposed at the interface with the grip plates, 
thus restraining all rigid body modes. 

The final item to be discussed in this section is the relationship of stresses, 
strains, and stress intensities in the remaining three- fourths of the specimen to those 
in the quarter model. In order to avoid confusion, assume that the upper quarter of 
the specimen, O^Q'^'TT/2, has been modeled and the two sets of boundary conditions, 
corresponding to the fundamental cases ,T a + b" and "a - b n , analyzed. Thus all answers 
will be given explicitly for the first quadrant. 

In the following, it is convenient to use a cartesian coordinate system as shown 
in Figure 8. Displacements in this system are related to those in the polar system via 
U r “ c °s 9 + U sin 0 

r X y (22) 

.U Q = JU y COS e - U x Sin 0 
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Figure A-8. - Cartesian coordinate system. 
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Applying these transformations to equations 13 shows that 


AJ (r, 0 +7T) = -JU { r, 6) 

X X 

JJL (r, 0 +7r) '= -U y (x> e) 


If one now introduces a second cartesian system that is rotated by'TT from the 
first then the displacements of the two points in their own systems are equal. 
Thus stresses and strains are also equal. Now stresses and strains are second 
order tensors and the transformation associated with the Tf rotation, or inversion 
of axes, is an identity transformation, thus stresses and strains are equal when 
expressed in either system. 


£ xx ( r , 0 +TT ) 

£ (r, 9 +TT ) 

^ yy 

6 xy (r, 0 +-,r ) 


(T (r, 0 + Tr > 

cr yy < r ’ 0 + ' T) 

CT xy < r > 6 +7r > 


II 

rTN 

(r, 0) 

XX 


/- 

= tl 

(r, 0) 

yy 


= c 

xy 

(r, 0) 


= U 

(r„ 0) 

XX 


= <r 

(r, 0) 

yy 


=<r 7 

(r, 0) 

xy 



(24) 


Since stresses are equal, stress intensities are necessarily equal also 

K.(r, 0 + TT) = K.(r, 0) (25) 

These results may be expressed with the point designations used earlier as 

{e (p 1 ) j = |e<p>) {£ (q 1 )) -{£<©} 

{cr (p 1 )] - {(r <p)j - {cr (q 1 )} ={<r «>} ««) 

| KCP 1 ) } - { K(P) ) -/ K(Q X ) ]- - { K(Q) } 

where j € (P) f denotes the column vector of strains evaluated at point "P", etc. 
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Thus the response in quadrants 2 and 3 is identical to that in 4 and 1 respectively 
for all loadings and they (2 and 3) need no further discussion* 

In considering quadrants 1 and 4, it is necessary to treat the two fundamental 
load cases separately. Consider the case T, a + b 11 first. Transforming equations 17 
via equations 22 shows that 

^ x ^ a+b = X ^ a+b 

(27) 

y (C ^a+b = ' W y (P) a+b 


Hence a second cartesian system obtained by inverting the "y" axis in the reterence 
system will render the points equal in their individual systems. The tensor trans- 
formation between this pair of coordinate systems produces a sign change on shearing 
terms, thus 1 




1 OD 

i 

i 

fe (?)") 

[ ^ XX j 

{ £ yy (Q) ' 

f _ 

0 lO 

i 

e yy ff> ■ 

^ £ xy^j 

a+b 

0 0-1 

1 

f (P) 
t' 1 ' xy v J 


(28) 


and the same relationship is applicable for stresses. Now and are related 

to normal and shearing stresses respectively, thus 

; 

f K l (P) 


( k i«>1 1 

" 1 o 

f 

I*** (Q) J a+b 1 

J 

1—1 

I 

0 

V. 


f 


K 


II 


(P)j 


(29) 


a+b 


For the case "a - b" transforming equations 21 gives. 


* x (Q) a-b = -^x <P >a-b 

(30) 

U (Q) , = h 

y a-b y a~b 
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hence inversion of the ri x ,! axis will make the points equal in their own system* 
The transformation between these two systems changes the sign of normal stresses, 
thus 




" -l oo~ 

1 


e „«> 

- 

0-10 

' 

€ < p ) 
yy 

l e >>j 

a-b ' 

o O 1 


£ (p) 

v x y J 


(31) 


and 


f K Z <Q» 

K IZ (Q)j 


a«b 


-1 O 

O 1 


\ (P)l 


K n ( p )i 


a-b 


(32) 


IV. TYPICAL FINITE ELEMENT MODELS 

Figures 9 and 10 show typical finite element models for the straight and branched 
crack configurations. In both specimen configurations, the region 0— 0 -Tf/2 was 
chosen for modeling. The center of the crack is located at the origin. In Figure 9 
the crack extends to the center of the semicircle. In Figure 10 it follows the "x" 
axis to the point shown and then follows a straight path to the center of the circle. 
The semicircular and circular regions are those occupied by the special crack tip 
elements. Half an element is used in the straight crack quarter model since the speci- 
men is' symmetric. Constant strain triangular elements were used in modeling the 
remaining region with small elements near the corner cuts for stress concentrations. 

A uniform displacement was applied at locations corresponding to the bolt holes 
to simulate rigid grip plates. Loads required to produce these displacements were 
determined from the analysis and used in normalizing results. Several different 
models were analyzed to investigate the effects of flow length and branching angle. 
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Figure A-9. - Typical finite-element model, straight crack. 
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Figure A-10. - Typical finite-element model* branched crack 





V. NUMERICAL RESULTS 


Normalized stress intensity factors and crack center displacements are shown 
in Figures 11-14. All calculations are for the case of plane strain (7 / = 1 / ) 
and V = .33 throughout. Normalized values were obtained by dividing the computed 
results by the corresponding value for an infinite plate in plane strain with a 
crack oriented at 45° from the applied loads and under the appropriate bi-axial 

i ■ / / 

loading (a + b or a - b) . For the branched crack cases the crack length used in 
normalizing was ’the sum of lengths as shown below. 



In referring to crack length, no dimensions are used. Any consistent dimensional 
system may be used throughout and the results will be applicable if the proportions 
of crack length to overall specimen geometry are maintained. For example, a crack 
length of two (2a - 2, a = 1) is applicable for a crack length equal to 1/10 of 
the grip plate width, which is shown to be twenty in Figure 9. The overall specimen 
geometry shorn in Figures 9 and 10 was maintained throughout and the crack geometry 
varied. 
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a ("a -fa") _ 
v("a + b") 

V = 0.33 
Plane Strain 


2 3 4 

LENGTH, a 


for straight cracks. 







L 


Stress intensity factors vs. flaw length for the straight crack are shown in 
Figure 11. For this configuration, the "a + b" load case produced pure K^. a n d 
crack opening. The "a - b" case gives pure and crack shearing. It should be 
noted that "a + b" is hydrostatic tension and "a - b" is pure shear in the 45 
orientation. This figure shows that is essentially constant over the range 
considered and that increases with increasing flaw length. Combining the 

scaling relationships discussed earlier with superposition and normalization leads 
to the following equations for stress intensities of straight cracks under any 
load combination. 


k.j. — % ( vj y ^ a ^ 


K-J.J = %(CT v “ O' H ) Vrr a K. 


H' 


II 


(33) . 


where K and are given by Figure 11 and -T R are average vertical and 

horizontal applied stresses associated with an uniform displacement of the grip 
plates. 

Corresponding crack displacements are shown in Figure 12. They are also 

normalized to plane strain infinite plate values with ~7 / = .33. The displacements 

shown are for a point at the center of the crack and on the upper face; r = 0 + 

-f ^ 

and 0=U/2 in Figure 2 or x = 0 and y = 0 in Figure 8. u and v are normalized 
shearing and opening displacements respectively. These curves show that the nor- 
malized displacements increase slightly with increasing crack length. Crack 
opening displacements (COD) and crack shearing displacements (CSD) are twice these 
values as shown in the sketch below. 
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From superposition, symmetry, scaling and normalization, it follows that 


cod = «r v +r H ) V 


CSD . mszh& (l r- v -«r H ) u 


(34) 


Similar results fox~ the branched crack are shown in Figures 13 and 14. 

These results were also normalized to infinite plate values with the total crack 
length (a = a^ + a^) and main crack orientation (45° from applied loads). A 
single main crack length was investigated (a^ = 1) with two branch lengths 
(a^ 53 .5, 1) and three angles (<>' = 36°, 54°, 72°). 

These curves show several interesting features. Both modes, and K^, 
are present in both the hydrostatic ("a 4- b !t ) and shear ( ?r a - b M ) loadings. 

Opening stress intensities are significantly larger than shearing intensities and 
the shear loading produces a large closing intensity. The branch length effect 
is nearly square root; as accounted for in the normalization. The effect of branch 
angle (oc ) is greater on than K^„, however the is smaller over the range 
studied. Straight crack results are also shorn along the ordinate for reference 
but these points have not been connected as the two problems are basically different. 
Consideration of superposition, symmetry; scaling and' normalization leads to equa- 
tions for stress intensities as in the straight crack case. However both branches 
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must be treated here as they have different intensities. Let the branches be 
denoted as shown below: t 


P 



1 1 

Intensities at P and Q will be equal to those at P and Q respectively, for any 
loading, as shown by equation 13. Intensities at P and Q are given by 


(k y 
1,2 

j 

* T V + ' r H 

! ; _ 

| ; V i 


f 

j 
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K 

K i,Q 

; = %\ ;, s7 a 

1 

'f v + j h 

| 

j ! - (r v -*v ; 

1 ! ^ 

K ’ 

1, 

i 

1 

i 

■ ( ' r V +<r H ) j j 

1 ! < 


i ) v “ J ^ 


V 



(35) 


were P and Q denote branches, "a + b f1 and "a - b TI denote the basic load cases. 
These equations show that the two branches have different intensities for any 
general load condition. 
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Normalized crack displacements for branched cracks are shown in Figure 14. 

These curves show that opening displacements (l/) are approximately equal to 
infinite plate values but shearing displacements are greater and they increase 
with increasing branch angle. Equations for COD an'd CSD are the same as for 
straight cracks, Equation 34, since the specimen is still symmetric. 

These four sets of curves, Figures 11-14, are the main results, from this 
analysis and they are applicable for any load combination, thickness, and Young's 

l 

modulus. These curves, along with Equations 33-35 permit calculation of stress 
intensity factors and crack displacements. In applying these results, the specimen 
configuration should be .as shown in Figures 9 and 10, and the specimen should be 
sufficiently thick to produce plane strain- with a material having a Poisson's 
ratio of .33. Several other results from the analysis are interesting and lead 
to a better understanding of the mixed mode problem. 

Figure 15 shows normalized crack displacements for a straight crack under the 
hydrostatic loading ("a + b"). Computed results were normalized by the infinite 
plate value at the crack center. The opening center value is slightly greater 
than an \nfinite plate and the distribution is ellipsoidal like an infinite plate. 
Sheari'ng displacements are very small over the entire crack length. Adjacent 
points on opposing crack faces have equal shearing displacements and equal but 
opposite opening displacements. Similar results for the shear loading ("a - b") 
are shown in Figure 16. These results are normalized by the infinite plate displace 
ment parallel to the crack and at its center. In this case, displacements perpen- 
dicular to the crack are of the same magnitude as those parallel to the crack. 
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Adjacent points have the same perpendicular displacement and equal but opposite 
parallel (shearing) displacements. The shearing displacement distribution is 
ellipsoidal as is an infinite plate. Displacements for other crack lengths x*ere 
similar to those shown here. 

Crack displacements for a branched crack are shown schematically in Figures 
17 and 18. The points labeled TI A, B, C M are at the branch juncture. These are 
distinct points in the finite element model but they all have the same coordinates. 
Points "D, E TT and "F, G ,f are similar and located where the crack intersects the 
special element. Crack tips are located at rT P, Q" and their displacements are not 
shown in these figures. Figure 17 shows that the ,f a + b fT loading produces opening 
of the main crack and both branches. Figure 18 shows sliding in the main crack, 
opening in the "Q" branch, and closing in the "P Tt branch. Closing is greatest 
between the points "A" and M C ,r . This closing is consistent with the negative 
shown in Figure 13. In the finite element model adjacent points on crack faces 
have the same coordinates and they are connected only via the elements; there 
is no direct connection across the crack faces. Thus the analysis permits them to 
overlap under load as shown on the TI P TT branch. This is physically impossible 
(impenetrability of matter) , however several other things must be considered before 
deducing that contact will occur in a test. First the relative porportions of 

n »y y* yj 

the two basic loadings (a + b and a - b) must be considered for the loading applied 

(CT_ and O' _) . If opening from the 11 a + b 11 component exceeds closing from the "a - b !I 

V H 

component then contact will not occur. Secondly, any residual opening such as 
from residual stresses, plasticity, finite width saw cuts, etc., must be considered 
as they were, not included in the analysis. 
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One model was reanalyzed under the assumption of plane stress to compare 
plane stress with pi are strain* It was a straight crack configuration with a 3,2 1.5 
and the stress intensities arer 


n 

a 4- b 

Plane Stress • 1*111 

Plane Strain 1.115 


a - b 

1.043 

1.088 


Both sets of stress intensities are normalized to infinite plate plane strain 
values, thus they are directly comparable. The differences are very small and 
probably less than the accuracy of the analysis. 
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APPENDIX B 

PRECRACKING STUDY AND PARAMETERS 
PRECRACKING STUDY 


Development trials were run to determine the relation between 
pressure level and crack growth for pressure-loaded flaws . Data 
from these trials are given in Table B— 1. The development specimens 
were 15% x 8.9 x 1.14 cm (6 x 3% x 0.450 in.) with the flaw parallel 
to the long direction. All specimens were cut from the same sheet 
of 2219-T87 material. A 0.635-cm (%-in.) hole was drilled through 
the specimen and used as a starter for a 5-cm (2— in.) long saber 
saw cut. Notches were put in the ends of the cut with a sharpened 
saber saw blade. The notches were scored with a razor blade, each 
corner of the blade being used only once. Cycling pressure was 
supplied by the hydraulic system of an MTS machine at 7 cycles/ 
second. Two of the best examples from the development series are 
shown in Figures B-l and B-2. The lightest areas are those of 
crack growth; the darkest are the original saw cut flaws and the 
saw cuts after cracking to permit breaking open of the specimens; 
the intermediate are the fracture areas. 


SPECIMEN PRECRACKING 

Based on the experience gained during the precracking develop- 
ment effort, initial values of pressure and number of cycles were 
established for the test specimens. However, in every instance of 
specimen precracking the precracker was opened several _times to 
permit tracking of the progress of the sharpening. After the first 
few specimens, the precracking effort became somewhat routine and 
complete data sets were not kept. The available data on pre- 
cracking the specimens is contained in Table B-2. 
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TABLE B-l, - PRECRACKING DEVELOPMENT TESTS 


Specimen 

designation 

Flaw orientation 
with respect to 
grain, deg 

Razor 

blade 

scoring a 

Maximum^ 

Number 

of 

cycles 

Flaw growth 

each end 

Flaw growth character 

pressure 

mm 

mils 


A 

0 

0 

B 

3 200 

0.25 - 1.91 

10 - 75 

Concentrated at one end. Uneven, out of 
plane. 

B 

0 

B 

A 

100 000 

0.38 - 1.91 

15 - 75 

Uneven at both enas, out of plane. 

C 

45 

B 

A 

61 000 

3.56 - 3.30 

140 - 130 

Uneven, out of plane. 

D 

45 

0 

B 

60 000 

3.81 - 0 

150 - 0 

In plane, one end only. 

E 

45 

0 

A 

76 000 

3.18 - 3.18 

125 - 125 

Out of plane. 

1 

45 

0 C 

B 

125 000 

1.65 - 1.65 

65 - 65 

Out of plane and uneven. 

2 

45 

0 

. B 

10 ooo 

0-0 

0-0 

No growth. 




A i 

52 000 

2.28 - 1.78 

90 - 70 

Out of plane and uneven. 

3 

45 

0 

l 

A i 

65 000 * 

3.30 - 2.54 

130 - 100 

In plane. 

4 

45 

0 d 

A 

40 000 

1.91 - 2.54 

75 - 100 

In plane one end, out of plane on other. 

1A 

45 

B d 

A 

45 000 

2.28 - 3.81 

90 - 150 

Out of plane and uneven. 

2A 

45 

B d 

"A 

25 000 

4.06 - 3.81 

160 - 150 

Out of plane. v 

3A 

45 

B d 

B 

20 000 

7.62 - 8.89 

300 - 350 

Even and in plane. 

4A 

45 

B d 

B 

10 000 

5.33 - 4.57 

210 - 180 

Even and in plane. 

5A 

45 

B d . 

C 

5 000 

4.32 - 4.57 

170 - 180 

Even and in plane. 


a O - One draw through each end* B - Two draws through each end, one from each side. 
b A - 2070 N/cm 2 (3000 pai); B - 2400 N/cm 2 (3500 psi)j C - 2750 N/cm 2 (4000 pai). 

C Exacto blade. 

Resharpened saber saw blade. 

Note ; All specimens: 2219-T87, t - 1,14 cm (0.450 in.), cut from one sheet. Initial flaw length - 5 cm (2 in,), flaw cut with 

saber saw. End V-notches made with sharpened saber saw. Notches scored with razor blade. Cycle rate ° 7 cycles/ second. 





2400 N/cm 2 (3500 psi), 20 000 Cycles. 



Figure B-2 . - Precracking Development Specimen 4A, 
2400 N/cm 2 (3500 psi), 10 000 Cycles. 
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TABLE B-2 


PRECRACKING PARAMETERS 


Specimen 

number 

Face 

Precracking flaw growth 

Average flaw 
length after 
sharpening 

Precracking parameters 

End A 

End B 

Cycles, 

thousands 

Press 

ure 

mm 

mils 

mm 

mils 

cm 

in. 

N/cm 2 

ksi 

1 

A 

1.7 

70 

1.1 

45 

6.541 

2.575 

5.5 

24.1 

3.5 


B 

1.0 

40 

1.7 

70 






2 

A 

1.7 

70 

1.2 

50 

5.385 

2.120 

13 

24.1 

3.5 


B 

2.0 

80 

1.2 

50 



1 

27.6 

4.0 

3 

A 

1.8 

75 

1.7 

70 

5.469 

2.153 

5 

27.6 

4.0 


B 

1.5 

60 

1.7 

70 



8 

24.1 

3.5 

4 

A 

1.8 

75 

1.6 

65 

5.423 

2.135 

5 

27.6 

4.0 


B 

1.6 

65 

1.6 

65 



4 

24.1 

3.5 

5 

A 

1.2 

50 

1.7 

70 

5.283 

2.080 

6 

22.1 

3.2 


B 

0.8 

30 

1.5 

60 






6 

A 

1.8 

75 

1.0 

40 

4.928 

1.940 

5 

22.1 

3.2 


B 

1.0 

40 

1.5 

60 






7 

A 

3.2 

130 

5.6 

230 

5.166 

2.034 

135 

22.1 

3.2 


B 

1.0 

40 

3.2 

130 



197.8 



8 

A 

1.8 

75 

1.5 

60 

6.375 

2.510 

112 

34.5 

5.0 


B 

3.7 

150 

3.2 

130 






9 

A 

2.5 

100 

3.4 

140 

5.474 

2.155 

50 

34.5 

5.0 


B 

4.4 

180 

6.6 

270 






11 

A 

1.7 

70 

1.6 

65 

10.60 

4.175 

8 

24.1 

3.5 


B 

1.8 

75 

0.5 

20 



20 

27.6 

4.0 

12 

A 

2.1 

85 

2.3 

95 

10.73 

4.225 

20 

27.6 

4.0 


B 

2.5 

100 

2.1 

85 






14 

A 

1.7 

70 

1.8 

75 

5.029 

1.980 

40 

34.5 

5.0 


B 

2.9 

120 

2.6 

105 






15 

A 

1.3 

55 

0.1 

5 

10.26 

4.040 

20 

17.2 

2.5 


B 

1.5 

60 

1.3 

55 



60 

24.8 

3.6 

16 

A 

4.2 

170 

3.9 

160 

10.92 

4.300 

25 

24.8 

3.6 


B 

4.9 

200 

4.9 

200 






17 

A 

2.0 

80 

2.7 

110 

10.49 

4.130 

14.6 

27.6 

4.0 


B 

1.8 

75 

1.3 

55 






18 

A 

1.2 

50 

0.6 

25 

10.34 

4.070 

6 

22.4 

3.3 


B 

0.4 

15 

0.5 

20 



9 

27.6 

4.0 

23 

A 

1.2 

50 

3.2 

130 

10.57 

4.160 

65 

22.4 

3.3 


B 

3.3 

135 

2.5 

100 






24 

A 

1.5 

60 

1.7 

70 

10.52 

4.140 

60 

22.4 

3.3 


B 

0.2 

10 

0.1 

5 






25 

A 

0.1 

5 

0.2 

10 

10.17 

4.005 

40 

22.4 

3.3 


B 

2.3 

95 

1.3 

55 






26 

A 

0.1 

5 

4.8 

195 

10.72 

4.220 

a 

24.1 

3.5 


B 

4.0 

165 

1.5 

60 






27 

A 

1.0 

40 

1.1 

45 

7.701 

3.032 

35 

24.1 

3.5 


B 

0.1 

5 

0.2 

10 






data. 
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APPENDIX C 


MIXED-MODE SYSTEM CHECKOUT SUMMARY REPORT 


Note ; This appendix is a copy of a report prepared by the Denver 
Division for Corporate Headquarters, Martin Marietta Corpo- 
ration. The checkout was performed as part of the facility 
procurement program of the Denver Division, and not as part 
of the mixed-mode program. Calibration data in Figure C-l 
are incorrect for the machine in the configuration ulti- 
mately put into service. 
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MIXED MODE SYSTEM - CHECKOUT SUMMARY REPORT 


System Operational Verification Objective : The objective of this phase of 

the checkout procedure was to verify the functionality of the Mixed Mode 
Test Machine and associated controls operating in concert with the Ur lb. 
MTS machine and to develop operational procedures. 

Results 


Note: Paragraph numbers used below refer to the corresponding objective 

section in the Checkout and Acceptance Plan. 

5.0 MMTM Uniaxial Checkout 

5.1.1 Horizontal axis load carrying capability was demonstrated at 300 kips 
in tension and 100 kips in compression. These loads exceed the load require- 
ments for contracted testing which can be forseen. Strain measurements 
obtained on the checkout specimen indicated that due to stress concentrations 
in the specimen corner the maximum tensile load which could be carried on the 
.3 inch thick 4340 steel specimen is approximately 140 kips rather than tl^e 
originally estimated 683 kips. The maximum compression load applied of 100 
kips likewise exceeds the compressive load requirements of the Mixed Mode 
Contract. 

5.1.2 The MMTM "load cell" was calibrated for loads from -100 kips to 
+300 kips. The resulting calibration data is shown in Figure 1. 

5.1.3 The strain distribution resulting from MMTM loading only is shown 
in Figure 3. Strain gage locations are shown in Figure 2. 

5.1.4 The dynamic response of the MMTM to cyclic ramp loads from function 
generator inputs was evaluated as a function of frequency at + 25, +75 
and + 100 kips. It was determined that the maximum frequencies which 
could be used were 1.0, 1.0, 0.9, and 0.5 Hz, respectively. The limiting 
factor in each case being an apparent resonance of the hydraulic lines. 

Checkout of the computer controlled operation of the MMTM was per_ 
formed using the programmed load input shown in 2.1 of the checkout " 
plan. The cyclic rate used was 8 sec/cycle. 

6.0 MMTM Support System - Final Check 

With the specimen pinned to both sides of the MMTM and the actuator 
side of the MTS and the MMTM supported on the three support jacks the 
MTS was cycled through + .125 inch stroke at frequencies up to 1 Hz.' The 
strain gages monitored during this test essentially followed the cyclic 
input (Figure 4) - indicating acceptability of the support system. 

Inertia effect of transverse fixture /during cycling of MTS machine intro- 
duces only minimal bending strains in the specimen. A 1/8 inch stroke 
at 1 cycle/ sec, which is larger and faster than will ever be used, pro- 
duced not over 700 psi stress at the most critical locations for detecting 
in-plane bending. 
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7*0 


MTS Checkout-Uniaxial 


7.5*1; 7.5.2 The strain distribution due to uniaxial loading in the MTS 
axis was determined incrementally at + 100 kips (Figure 5). 

7.6 Cyclic functioning of the MTS was verified using function generator 
inputs . 

7.7 Operation of the MTS system using computer control was verified 
using the load history shown in 2.1 of the checkout plan. 

8.0 Biaxial System Checkout 


Prior to proceeding with the programmed biaxial loading static biaxial 
loads were applied to evaluate strain distribution. The results of 
loading to + 100 kips in 25 kip increments is shown in Figures 6a and> b* 

The strain distribution in biaxial tension is uniform to less than +5% 
along both axes up to 100 kips for a 9 inch diameter around the center of 
the specimen. In compression the strain distribution is uniform to with- 
in +5% at loads up to 75 kip for a -9 inch diameter and within a inch 
diameter at 100 kip. This uniformity of strain is well within the 
acceptance limits. 

The system was operated in the biaxial mode using computer control. 

The programmed load input was as shown in 2.1 of the checkout plan. Response 
of the system operating at 8 sec/cycle is shown in Figures 7a and 7b. In 
both figures the feedback from the load cells on each axis is shown. Both' 
axes exhibit some looseness going through zero. At this cyclic rate it is 
acceptable since the peak load levels achieved are not affected. The 
figures also show the response of strain gages located inside a 4% inch 
radius circle near the center of the specimen. Figure 7a shows the 
response of gages aligned with the MTS axis; Figure 7b those aligned with 
the MMTM axis. Two points with respect to the strain gages response should 
be noted: 1) In both cases the fidelity of the gage response reflects 

that of the load input and 2) The strains peaks agree with those pre- 
dicted analytically. 

The system was also operated at a cyclic rate of 4 sec/cycle 
(Figure 8). At this cyclic rate the feedback from the MMTM axis 
was distorted and load and strain peaks affected by reverberations in 
the system. It has been concluded that 4 sec/cycle is too fast to 
obtain meaningful test data. 

Conclusions 


1. The load carrying capability of the MMTM has been demonstrated to 
300 kip which is more than adequate for the forseeable test usage. 

2. The stress distribution in the specimen is extremely uniform under 
equal biaxial loads. 

3* System functions as required under computer control input to both axes. 
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on Opposite Side. 

20 Channels Total. Not to Scale. 


ONE INCH THICK ALUMINUM CHECKOUT SPECIMEN 

Figure 2 
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APPENDIX D 


TABLES OF PRIMARY CYCLIC AND STATIC DATA 


This appendix contains a description of the experimental re- 
sponse of each specimen and tables of the primary load and crack 
length data recorded during each test in, the program. All lengths 
in the tables are actual crack lengths. Where projected crack 
lengths were required in the analyses, the crack lengths in the 
tables were transformed to the appropriate projected length. 

SN-1 . - This sample was subjected to approximately 1750 cycles 
of loading in the pure K^. mode, load type A (equal, synchronized 

loads on each axis). Load was applied at very low stress intensity 

values | less than 5500 kN cm 2 /cm 2 or 5000 psi-in. 2 ) to check the 
operation of the recording equipment. Cycle testing of the speci- 
men then began at a stress intensity level of approximately 

11 000 kN cm 2 /cm 2 (lO 000 psi-in." 2 ), which corresponded to a maxi- 
mum biaxial load of 289 kN (65 kips). The table for SN-1 in this 
appendix lists the subsequent loads, the number of cycles at each 
load level, and the resulting crack length. The crack did not 
grow straight. However, each increment of crack growth was meas- 
ured as a straight line from the last position of the crack tip at 
which a measurement was taken. The values of crack length in the 
table are the average of the four values of crack growth measured 
at each end of the crack on each face of the specimen, added cumu- 
latively to the original flaw length. 

The 2673 cycles at 311-kN (70-kip) maximum load were applied 
during an effort to locate a noise on the horizontal axis that 
was causing a bump in the horizontal load feedback signal. The 
noise was localized to the bottom actuator on the transverse axis. 
Jamb nuts were installed on the actuator piston rod but the prob- 
lem persisted throughout the program. 

In the table, the 1102 cycles of loading at 405 kN (91 kips), 
primarily on the horizontal axis, was an attempt to cause the 
crack to turn back toward the 45-deg line. It appeared from the 
intersection of the crack leading edges with the specimen sur- 
faces, that the attempt was successful. Subsequently, the Tt phase 
shifting" technique previously described was employed. However, 
the crack resumed its growth in the horizontal direction. 

The static test performed after completion of cycling did not 
yield usable data. The crack under the top loading tab in the photo- 
graphs in Appendix E was caused during specimen removal. 
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SN-2 . ' - This sample was cyclically tested under load type A 
as a duplicate of SN-1. The maximum load imposed during the 
cyclic tests was 592 kN (133 kips). After cyclic testing, SN-2 
was tested statically to fracture under equal biaxial loads . 
High-speed movies revealed that subcritical crack growth started 
•at the upper end of the crack at approximately 594 kN (134 kips) 
and at the lower end at approximately 604 kN (136 kips) . Frac- 
ture occurred at approximately 774 kN (174 kips) after subcritical 
growth of between 0.76 and 2.54 cm (0.3 and 1.0 in.), depending 
on which end of the flaw and which face of the specimen were con- 
sidered. Turning the crack appears to have started at the very 
beginning of cyclic testing and continued until the cracks were 
parallel to the final rolling direction. This is -the same 
phenomenon as observed in the test of SN-1. 

SN-3. - This specimen was the first cyclically tested under 
the mixed-mode condition (load type B) . Branching was immediately 
exhibited at both ends of the flaw and on both faces of the speci- 
men. All legs of the flaw continued to grow until the test was 
stopped. The test was stopped during cyclic testing after approxi- 
mately 22 700 cycles when there was sudden growth (3.8 to 5 cm or 
1% to 2 in.) of the horizontal legs of the cracks during one cycle 
of loading at the 489-kN (110-kip) load level. At the time this 
appeared to be a considerably premature failure. There was no 
further testing, cyclic or static. 

Sn- 4. - This specimen was cyclically tested under load type 
C for approximately 32 000 cycles up to a maximum load of 569 kN 
(128 kips) then tested statically to failure. Failure occurred 
under static loading at 850 kN (191 kips) as evidenced by sudden 
growth of the cracks into the grips. Branching of the flaw ap- 
peared to start immediately on cyclic testing. There was con- 
siderable subcritical crack growth of the horizontal legs of the 
branched cracks during static testing. The static load was ap- 
plied normal to the rolling direction. Due to operator error the 
load was not removed after failure but allowed to increase. At 
889.6 kN (200 kips)- the top loading tab of the specimen broke 
off It was thus learned that the central portion of the specimens 

must be weakened by sufficient flaw growth (or longer initial 
flaws for the specimens to be tested only statically) to the point 
where fracture can occur there and not across the tab notches . 


SN-5. _ This sample was cyclically tested under load type B 
for approximately 43 100 cycles up to maximum loads of 265 kN 
(59.5 kips) on one axis and 113 kN (25.5 kips) on the other, then 
tested statically to failure. Failure occurred at 545 kN (122.6 
kips) on the vertical axis with 230 kN (51.8 kips) on the hori- 
zontal axis. Branching of the flaw appeared to start immediately 
on cyclic testing. The branches parallel to the rolling direc 
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tion grew much faster than those perpendicular to it during the 
cyclic test. 

SN-6 . - This specimen was cyclically tested under equal bi- 
axial loads, load type A, for approximately 55 000 cycles up to 
a maximum load of 244 kN (54.8 kips). It was then tested stat- 
ically to failure under equal biaxial loads. Failure occurred 
at about 513 kN (115.4 kips). Crack growth during cyclic test- 
ing was straight until a shear lip started forming on the lower 
end of the flaw. On one surface of the specimen , the crack re- 
mainted straight while on the other it curved sharply to the 
horizontal. At fracture the upper end of the flaw grew vertically 
while the lower end grew horizontally. 

SN-7 . - SN-7 was cyclically tested under equal biaxial loads 
for approximately 40 500 cycles up to a maximum load of 512 kN 
(115 kips). Fracture occurred at the 696th cycle at 512 kN (115 
kips) . The length of the crack at fracture is unknown because no 
movies were taken during cyclic testing. It was decided to test SN- 
7 cyclically to failure to obtain as complete a plot of AG versus 
Aa/AN as possible. The flaw turned parallel to the rolling direc- 
tion soon after the beginning of cyclic testing and grew into the 
grips at fracture. The length of the crack at fracture has been 
estimated by extrapolating the compliance plot to the value of com- 
pliance measured on the last cycle, which was recorded on the oscil- 
lograph- The last cycle has been considered a static test in the 
data analysis. 

SN-8 . - SN-8 was cyclically tested under load type B for 
23 673 cycles at which time one of the load tabs broke off. The 
maximum cyclic loads achieved were approximately 514 kN (115.5 
kips) on one axis with 218 kN (49 kips) on the other. Though 
branching of the crack was apparent at the very beginning of 
cyclic loading, the vertical legs of the branches developed only 
slightly. Most crack growth took place in the horizontal direc- 
tion 5 parallel to the rolling direction. Vertical crack growth 
has been ignored in the data analysis as insignificant. 

SN-9 » - This was essentially a retest of SN-8 except that the 
final rolling direction was oriented vertically instead of hori- 
zontally in the testing machine. SN-9 was loaded cyclically for 
approximately 39 306 cycles. It was taken to failure under cyclic 
loading to provide a good comparison with specimen SN-7. The max- 
imum loads just prior to failure were approximately 503 kN (113 
kips) on one axis with 224 kN (50 kips) on the other. The crack 
branched and, in the early portions of the test, the branches 
grew at about the same rate. However, the vertical branches 
started to grow faster about half way through the test and growth 
of the horizontal branches almost ceased. Fracture was charac- 
terized by sudden growth of the vertical branches into the grips. 


171 



The length of the crack at fracture has been estimated by extra- 
polating the crack growth rate from the oscillograph records* The 
last cycle has been considered a static test* 

SN-ll. - This specimen was tested Statically with a compres- 
sion load on the vertical axis equal to about 40% of the tension 
load on the horizontal axis (load type D) . The initial flaw was 
approximately 10— cm (4— in.) long. Fracture occurred with 256— kN 
(57.6-kip) compression on the vertical axis and 615-kN (138.2- • 

kip) tension on the horizontal axis = 2.6j. There was 

some subcritical crack growth, starting at —153 kN (—34.4 kips) 
and 371 kN (83.4 kips) on the vertical and horizontal axes, . 
respectively. Crack growth was straight from the ends of the 
flaw directly into the comers of the specimen. There was no 
shear lip or turning of the crack. In the photographs of 
Appendix E, note the permanent shear deformation at the hole 
used for mounting the biaxial COD gage holder. 

SN-12 . - SN-12 was tested statically with compression on the 
vertical axis numerically equal to the tension on the horizontal 
axis, resulting in pure shear on the 10— cm (4— in.) 45— deg flaw 
(load type E) . Fracture occurred when the specimen buckled 
under the vertical compressive load in spite of the presence of 
the antibuckling plates. Simultaneously one of the loading tabs 
on the horizontal axis broke off. The fracture was vertical, 
perpendicular to the rolling direction, but also perpendicular 
to the tension load. 

SN-14. - SN-14 was cyclically tested under load type B; 
however, the shear was not reversed on each cycle. Approximately 
39 300 cycles were applied with the higher load, 245 kN (55 kips), 
always in the vertical direction, resulting in horizontal crack 
growth of about 0.76 cm (0.3 in.) off each end of the original 
flaw. Then the shear was reversed by putting the higher load 
of 245 kN (55 kips) on the horizontal axis while the vertical 
load was reduced to 89 kN (20 kips) . No crack growth was noticed 
for about 60 500 cycles. The horizontal load was increased to 
311 kN (70 kips) with the vertical load at 116 kN (26 kips) and 
vertical branches started to grow. On one side of the plate they 
grew from the ends of the original flaw. On the other side the 
top one grew from the original flaw but the bottom one grew from 
the end of the horizontal branch. Total cycles at the 311-kN 
(70-kip) level were approximately 24 500. The vertical branches 
had grown about 0.76 cm .(0.3 in.). The test was terminated and 
SN-14 was subsequently tested statically. It is interesting to 
note that once the horizontal branches were started, applying the 
high horizontal load did not induce vertical branches as it would 
have had the high load been alternated between -the vertical and 
horizontal axes on each cycle. SN-14 was tested statically under 
a vertical tension load. Fracture occurred at 1468 kN (330 kips) 
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with the horizontal branches running toward the grips then -turn- 
ing parallel to the rolling direction. The fracture surface ex- 
hibited large shear lips. 

SN-15* - This specimen was tested statically under equal 
biaxial loading |pure The initial flaw was approximately 

lO-cm (4-in.) long. Fracture occurred at 965 to 968 kN (217.0 
to 217.5 kips) on each axis after approximately 1.2 cm (0*5 in.) 
of sub critical crack growth, at each end of the flaw. The crack 
turned toward the rolling direction but was still generally at 
45 deg to the loading axis when the critical length was reached* 

SN-16 . - SN-16 was tested statically with loading only on the 
horizontal axis = 1 Again, the initial flaw was approx- 

imately 10 cm (4 in.) long. Fracture occurred at 1108 kN (249' kips). 
The crack propagated normal to the tension load, and therefore nor- 
mal to the rolling direction, but turned parallel to the rolling 
direction before entering the grips. 

SN-17 * - SN-17 was tested statically with a compression load 
on the vertical axis equal to about 40% of the tension load on the 
horizontal axis (load type D) . The initial flaw was approximately 
10-cm (4-in.) long. Fracture occurred with a 1486-kN (334-kip) 
tension on the horizontal axis and a 614-kN (138-kip) compression 
on the vertical axis = 2.5 j. There was some subcritical 

crack growth. Crack growth was straight from the ends of the flaw 
directly into the corners of the specimen. There was no shear 
lip or turning of the crack. In the photographs (Appendix E) , note 
the permanent shear displacement at the center hole. 

SN-13 . - SN-18 was tested statically with compression on the 
vertical axis numerically equal to the tension on the horizontal 
axis, resulting in pure shear on the 10-cm (4-in.) 45-deg flaw.- 
Fracture occurred with P Tr = -1103 kN (-248 kips) and P - 1130 kN 

"•V H 

(254 kips). There was some subcritical growth. Crack growth was 
straight into the corner notches. There were no shear lips. 

SN-23 * - SN-23 was tested statically with tension on the 
vertical axis and no load on the horizontal axis (load type C) . 

If there had been no subcritical growth, this loading arrange- 
ment would have produced a ratio - 1. However, the small 

amount of subcritical growth in the horizontal direction dras- 
tically changed the ratio at failure. Fracture occurred 
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at = 752 kN (169 kips) . Subcritical growth had started at 

690 kN (155 kips). In the photographs in Appendix E, the crack 
emanating from the lower left notch opened after the fracture. 

Since the fracture cracks ran into the grips, the specimen-grip 
assembly was still in one piece and vertical load was transmitted 
through the horizontal grip plates. Because the MTS testing 
machine was operated as a load-controlled device, it continued to 
apply load until the second crack resulted in excessive actuator 
travel. Specimen SN-23 should be compared with specimen SN-16. 

The only difference was that SN-23 had the load applied perpen- 
dicular to the rolling direction while the load was applied paral- 
lel to the rolling direction on SN-16. 

SN-24 and SN-25 . - These specimens were both tested statically 
with compression on the vertical axis numerically equal to the 
tension on the horizontal axis , resulting in pure shear on the 
flaws. The only -difference between the tests was that the tension 
load was applied perpendicular to the rolling direction on SN-25 
and parallel to the rolling direction on SN-24. However, there was 
little difference in the results. SN-24 failed at approximately 
1027 kN (231 kips) and SN-25 at 1005 kN (226 kips). There was 
some subcritical growth. Crack growth in both tests was straight 
into the corners. 

SN-26. - SN-26 was tested statically with a compression load 
on the vertical axis equal to about 40% of the tensxon load on the 
horizontal axis (load type D) . The initial flaw was approximately 
10-cm (4-in.) long. Fracture occurred with a 1281-kN (288-kip) 
tension on the horizontal axis and a 498-kN (112-kip) compression 
on the vertical axis 1T J Kj. = 2.4 ) when a loading tab on the 

horizontal axis broke off. There was some subcritical crack 
growth. Unfortunately, there is no load record to correlate with 
the small amount of crack growth that did take place. The tab 
failure initiated at a sharp flaw in one corner of the specimen 
that had gone unnoticed before the test. The data analvsis pre- 
sented earlier used the final crack length and maximum loads to 
compute a minimum G . 

Clf 

SN-27. - SN-27 was tested cyclically under a compression 
load on the vertical axis equal to the tension load on the hori- 
zontal axis on every cycle (load type E) . This loading was in- 
tended to produce a case of pure without shear reversal. 

However, the crack immediately .started to grow perpendicular to 
the tension load rather than in line with the 7.5-cm (3-in.) flaw, 
which meant that there was little shear on the crack tip but 
much tension. The cyclic testing was continued for a little more 
than 26,000 cycles until the crack grew into the grip areas. No 
static test was made since load would have been transferred 
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through the grip plates rather than through the fully cracked 
specimen. . In the photographs in Appendix E, the slant crack 
into the lower right corner was inadvertently caused when the 
specimen was being removed from the testing machine. 


STATIC TEST CRACK GROWTH FOR SN-6 


Vertical load, 

p v 

Horizontal load, 
P H 

Crack length, 
2a 

kN 

kip 

kN 

kip 

cm 

in* 

297 

66.8 

307 

69 

6.15 

2.42 

495 

111.2 

508 

114.2 

7.67 

3.02 

499 

112.2 

512 

115 

7.98 

3.14 

504 

113.4 

513 

115.4 

8.26 

3.25 

504 

113.4 

515 

115.8 

8.46 

. 3.33 

506 | 

113.8 

516 

116 

8.69 

3.42 

508 

114.2 

517 

116.2 

8.99 

3.54 


STATIC TEST CRACK GROWTH FOR SN-11 


Vertical load, 

p v 

Horizontal load, 
P H 

Crack length, 
2a 

kN 

kip 

kN 

kip 

cm 

in. 

-153 

-34.4 

371 

83.4 

10.61 

4.18 

-236 

-53 

569 

128 

11.20 

4.41 

-246 

-55.4 

594 

133.6 

11.99 

4.72 

-252 

-56.6 

606' 

136.2 

12.83 

5.05 

-•254 

-57 

612 

137.6 

13.28 

5.23 

-255 

-57.4 

614 . 

138 

13.64 

5.37 

-255 

-57.4 

615 

138.2 

13.89 

5.47 

-252 

-56.7 

615 

138.2 

14.27 

5.62 

-252 

-56.7 

616 

138.4 

15.90 

6.26 

-252 

-56.7 

615 

138.2 

16.41 

6.46 

-252 

-56.7 

615 

138.2 

- 17.17 

6.76 





Vertical load* 

V 

fcN 

kip 

323 

72.7 

399 

89.8 

426 

95.8 

466 

104.8 

491 

110.4 

512 

115.2 

527 

118.4 

535 

120.2 

539 

121.2 

542 

121.8 

544 

122.2 

544 

122.2 

545 1 

122.6 

545 

122.6 

545 

122.6 

545 

122.6 


Horizontal load , ’ 
P H 

Crack length, 
2a 

. . kN 



in. 

138 

31.0 

9.37 

3.69 

172 

38.6 

9.63 

3.79 

188 

42.2 

9.96 

3.92 

197 : 

44.2 

10.46 

-4.12 

207 

46.6 

10.62 

4.18 

215 

48.4 

11;20 

4.41 

222 

49.8 

11.63 ' 

4.58' 

225 

50.6 

11.91 

4.69 

227 

51 

12.12 

4.77; 

229 

51.4 

12.32 

’4.85' 

229 

51.4 

12.70 . 

5.00' 

229 

,51.4 

12.83 

5.05 

230 

51.6 

13.18 

5.19 

230 

51.6 

13.89 

5.47 

230 

51.6 

14.33 

5.64. 

230 • 

51.8 

15.52 

6.11 














STATIC TEST CRACK GROWTH FOR SN-4 


Vertical load, 

p v 

Crack length, 
2a 

kN 

kip 

cm 

in. 



10.74 

4.23 



750 

168.5 

10.92 

4.30 

772 

17.3.5 

11.05 

4.35 

794 

178.5 

11.13 

4.38 

816 

183.5 

11.38 

4.48 

841 

189 

11.84 

4.66 

843 

189.5 

11.91 

4.69 

845 

190 

12.12 

4.77 

847 

190.5 

12.57 

4.95 

850 

191 

13.34 

5.25 

850 

191 

13.67 

5.38 

850 

191 

14.10 

5.55 

850 

191 

14.38 

5.66 

Note: 

p H = 0. 




STATIC TEST CRACK GROWTH FOR SN-14 


Vertical load, 

p v 

Crack length, 
2a 

kN 

kip 

cm 1 

in. 

— 


7.52 

2.96 

1352 

304 

9.09 

3.58 

1443 

324.5 

9.50 

3.74 

1457 

327.5 

10.01 

3.94 

1457 

327.5 

10.29 

4.05 

' 1468 

330 

10.57 

4.16 

-P> 

CO 

330 

10.62 

4.18 

1468 

330 

11.56 

4.55 

Note’: 

f h- 




STATIC TEST CRACK GROWTH FOR SN-15 


Vertical load. 

Horizontal load. 

Crack length. 

r 

V 


P H 

2a 



kip 


kip 

cm 






10.26 

mm 







894 

201 

11.02 

Hi 


209 

930 

209 

11.30 

as 


214 

956 

215 

11.46 

4.51 

959 

215.5 

963 

216.5 

11.58 

4.56 

963 

216.5 

965 

217 


4.61 

965 

217 

965 

217 


-4.65 

965 

217 

967 

217.5 


4.75 

965 

217 

967 

217.5 


4.99 


177 


















STATIC TEST CRACK GROWTH FOR SN-16 















STATIC TEST CRACK GROWTH FOR SN-17 


Vertical load, I Horizontal load 
P„ 


Crack length, 
2a 


-558 

-125.5 

-574 

-129 

-583 

-131 

-594 

-133.5 

-592 

-133 

-601 

-135 

-601 

-135 

-603 

-135.5 


1350 303.5 

1390 312.5- 

1421 319.5 

1437 323 

1443 324.5 

1446 325 

1448 325.5 

1448 325.5 


10.49 4.13 
11.23. 4.42 
11 -.46 4.51 
11.84 4.66 
12.27 4.83 
12.47 4.91 
12.67 4.99 
12.80 5.04 
13.18 5.19 


STATIC TEST CRACK GROWTH FOR SN-2 


. load, 

r 

Horizontal load* 
P H 

Crack ; 
2 j 

kip 

kN 

kip 

cm 

133.5 

603 

135.5 

‘ 12.55 

152 

688 

154.6 

12.73 

156.6 

707 

159 

12.85 

159.4 

719 

161.6 

12.98 

163.2 

737 

165.6 

13.26 

166 

747 

168 

13.36 

167 

754 

169.4 

13.44 

1 7n A 


1 79 O 


171.4 

/ uu 

773 

173.8 

IJ * U_L 

13.84 

172.6 

776 

174.5 

14.05 

173 

778 

174.8 

14.30 

173.6 

778 

174.8 

14.55 

173.6 

781 

175.6 

15.54 




























STATIC TEST CRACK GROWTH FOR SN-18 


Vertical load, 

. p v 

Horizontal load, 
P H 

kN 

’kip . 

kN 

kip 


-676 -152 ' 681 153 

-1025 -230.5 1032 232 

-1079' ' -242.5 . 1083 '243.5 

-1103 -248 1101 247.5 

-1117 -251 1114 250.5 

-1123 -252.5 ' 1117 251 

-1125 -253 : 1119 251.5 

-1125 -253 ' .1121 '252 

-1130 -254 1121 252 

-1130 -254 1117 251 

-1130 -254 1108 249 


Crack length, 
2a 


11.91 
12.70 
13.56 
14.35 
14.48 
14.76 
15.44 

16.92 
17.75 


STATIC TEST CRACK GROWTH FOR SN-23 


Note: P„ = 0. 

H 


Crack length. 

2a 


cm 

in. , 

10.57 

4.16 

10.74 

4.23 

11.91 

4.69 


STATIC TEST CRACK GROWTH FOR SN-26 


Vertical load. Horizontal load. Crack length, 
P V P H 2a 


cm in. 

-112 | 1281 I 288 p 11.46 # 4.51 


Note 


loading tab broke off 


kip 

kN 

kips 

-112 

1281 

288 . 


L8 









































STATIC TEST CRACK GROWTH FOR SN-24 


Vertical load, 

p v 

Horizontal load, 
P H 

Crack length, 
2a 

kN 

kip 

kN 

kip 

cm 

in. 

-692 

-155.5 

692 

155.5 

10.52 

4.14 

-990 

-222.5 

994 

223.5 

11.40 

4.49 

-1010 

-227 

1014 

228 

12.04 

4.74 

-1028 

-231 

1023 

230 

12.67 

4.99 

-1028 

-231 

1025 

230.5 

12.90 

5.08 

-1032 

-232 

1025 

230.5 

13.28 

5.23 

-1034 

-232.5 

1025 

230.5 

13.56 

5.34 

-1034 

-232.5 

1023 

230 

13.82 

5.44 

-1036 

-233 

1023 

230 

14.50 

5.71 

-1036 

-233 

1021 

229.5 

15.47 

6.09 

-1036 

-233 

1008 

226.5 

16.48 

6.49 

-1036 

-233 

1008 

226.5 

18.01 

7.09 


STATIC TEST CRACK GROWTH FOR SN-25 


Vertical load, 

p v 

Horizontal load, 
P H 

Crack length, 
2a 

kN. 

kip 

kN 

kip 

cm 

in. 

-665 

-149.5 

665 

149.5 

10.19 

4.01 

-974 

-219 

976 

219.5 

11.38 

4.48 

-994 

-223.5 

992 

223 

11.84 

4.66 

-1001 

-225 

999 

224.5 

12.24 

4.82 

-1005 

-226 

1001 

225 

12.78 

5.03 

-1008 

-226.5 

1001 

225 

13.26 

5.22 

-1008 

-226.5 

999 

224.5 

13.59 

5.35 

-1010 

-227 

999 

224.5 

14.40 

5.67 

-1010 

-227 

996 

224 

16.64 

6.55 


-227 

974 

219 

17.45 

6.87 
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LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-1 


Nominal 

P , kN 

■max 

(kip.) 

Horizontal load, P^, 
kN (kip) 

Vertical load, P^., 
kN (kip) 

Increment 
of cycles , 
AN 

Crack length, 
2a 



Max 

Min 

Max 

Min 

cm 

in. 







2a. = 6.541 
i 

2.575 

<156 


* 



: 

1 740 

6.566 

2.585 

(<35) 




► 

* 



289 

294-300 

18-27 . • • 

289-294 

2'9-33 ■ 

7 756 

\ 6.777 

2. ’66 8 

C65> 

(66.0-67.5) 

(4. 0-6.0) • 

(64.5-66.0) 

(6.5-7. 5) 

8 737 

7.127 

2.806 

* 


. * 

* 

• 

8 255 

7.462 

2.938 

347 

348-360 

29-44 ' • 

347-351 

33-40 

2 481 

7.653 

3.013 . 

(78) 

(77.5-81.0) 

(6.5-10.0) • 

(78.0-79.0) 

(7. 5-9.0) 

3 620 

8.087 

3.184 






5 383 

8.588 

3.-381 

405 

409-414 

33 

400 

44 

2 417 

9.129 

3.594 

(91) 

(92.0-93.0) 

(7.5) 

(90.0) • 

(10.0) 

2 097 

9.563 

3.765 


* 

- 



882 • 


.... . 

311 

302 

22 

307 

40 

' 2 673 



(70) 

(68.0). 

(5.0) 

(69.0) 

(9.0) • 

• 



408 

'403 . 

56 

47 ’ 

47 

1 102 

9.959 

3;921 

(91) 

(90.5) 

(12.5) 

(10.5) 

(10.5) 




408 

408 

42 

403 ‘ 

42 

2 395 

10.44 

4.110 

. (91) 

(91.0) 

(9.5) 

(90.5) 

(9.5) 




463 

466-472 

41-48 

460-463 

42-45 

1 059 

10.90 

4.291 

(104) 

(104.8-106.1) 

(9.2-10.8) 

(103.3-104.0) 

(9.5-10.1) 

500 

11.27 

4.435 






404 

11.59 

4.564 

534 

540-542 

43 

530-536 

39-40 

209 

12.10 

4.765 

(120) 

(121.4-121.8) 

(9.6) 

(119.1-120.4) 

(8. 8-9.1) 

146 

12.46 

4.904 






203 

13.01 

5.121 

565 

576 

45-47 

560-562 

36-39 

105 

13.65 

5.373 

(127) 

(129.4-129.5) 

(10.1-10.5) 

(126.0-126.3) 

(8. 2-8. 8) 

49 

13.97 

5.499 






50 

14.21 

'5.596 





Total 

52 263 





















LOADING SCHEDULING AND CRACK GROWTH FOR CYCLIC TEST OF SN-2 


Nominal 

Horizontal 

load, P„, 

Vertical load, P , 

Increment 

Crack length, 

P M * 

Max 

n 

kN (kip) 

V 

kN (kip) 

of cycles, 
AN 

2a 


(kip) 

Max 

Min 

Max 

Min 


cm 

in. 







2a.. = 5.385 
i 

2.120 

222 

227-234 

20-27 

220-225 

24-27 

6 561 

5.458 

2.149 

(50) 

(51.0-52.5) 

(4. 5-6.0) 

(49.5-50.5) 

(5. 5-6.0) 

11 000 

5.644 

2.222 






12 030 

5.928 

2.334 

289 

282-285 

27-31 

287-285 

29-31 

5 022 

6.160 

2.425 

(65) 

(63.5-64.0) 

6. 0-7.0) 

(64.5-64.0) 

(6. 5-7.0) 

5 991 

6.449 

2.539 






4 969 

6.690 

2.634 

356 

376-378 

36-38 

351-356 

38-42 

3 056 

7.036 

2.770 

(80) 

(84.5-85.0) 

(8. 0-8. 5) 

(79.0-80.0) 

(8. 5-9. 5) 

3 015 

7.407 

2.916 






2 000 

7.689 


423 

416-425 

40-42 

416-420 

47-49 

1 020 

8.029 

3.161 

(95) 

(93.5-95.5) 

(9. 0-9. 5) 

(93.5-94.5) 

(10.5-11.0) 

1 000 

8.252 

3.249 






997 

8.547 

3.365 

489 

485-494 

49-58 

489-507 

49-62 

518 

8.791 

3.461 

(no) 

(109.0-111.0) 

(11.0-13.0) 

(110.0-114.0) 

(11.0-14.0) 

500 

9.050 

3,563 






500 

9.370 

3.689 






500 

9.733 ' 







500 

10.18 

4.008 

556 

578-592 

62-76 

578-583 

53-62 

110 

10.70 

4.213 

(125) 

(130.0-133.0) 

(14.0-17.0) 

(130.0-131.6) 

(12.0-14.0) 

ioo 

11.06 

4.354 






50 

11.39 

4.484 






100 

11.77 

4.635 






75 

12.12 

4.770 






75 

12.55 

4.940 





Total 

59 689 




K 1 


00 

LO 















LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF .SN-3 


Nominal 

Max’ 

(kip) 

Horizontal load, P_ T , kN (kip) 
rt 

Vertical load, P^, kN (kip) 

Increment 
of cycles, 
AN 

Crack length, 2a 

Vertical 

crack 

Horizontal 

crack 

High max 

Low max 

Min 

High max 

Low max 

Min 









cm 

in. 









2a ± - 5.469 

2a. = 2.153 

X 

222 

229-231 

85-89 

11-13 

220-222 

82-85 

9 

13 657 

5.735 

5.801 

(50) 

(51.5-52.0) 

(19.0-20.0) 

(2. 5-3.0) 

(49.5-50.0) 

(18.5 r 19.0) 

(2.0) 


(2.258) 

(2.284) 

311 

320 

100-122 

11 

309-311 

96-113 

9-11 

5 014 

6.027 

6.017 

(70) 

(72.0) 

(22.5-27.5) 

(2.5) 

(69,5-70.0) 

(21.5-25.5) 

(2. 0-2. 5) 


(2.373) 

(2.369) 

400 

407 

129-131 

9-11 

391-394 

100-122 

11-13 

2 740 

6.281 

6.337 

(90) 

(91.5) 

(29.0-29.5) 

(2. 0-2. 5) 

(88.0-88,5) 

(22.5-27.5) 

(2. 5-3.0) 


(2.473) 

(2.495) 

489 

489 

191 

22 

489 

187 

13 

1 243 

6.541 

6.546 

(110) 

(110.0) 

(43.0) 

(5.0) 

(110.0) 

(42.0) 

(3.0) 


(2.575) 

(2.577) 







Total 

22 654 






LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-4 


Horizontal load, 



Vertical load, 
P v , kN (kip) 


Increment 
of cycles, 
AN 


Crack length, 2a 


Vertical 

crack 


2a. = 2a. = 

i x 

5.423 2.135 


Horizontal 

crack 


231-236 

(52.0-53.0) 

276 

(62.0) 

316-325 

(71.0-73.0) 

374 

(84.0) 

400 

(90.0) 

449 

( 101 . 0 ) 

476 

(107.0) 

529 

(119.0) 

565 

(127.0) 


13-18 
(3. 0-4.0) 

13 

(3.0) 

13 

(3.0) 

18 

(4.0) 

18 

(4.0) 

18 

(4.0) 

22 

(5.0) 

13 

(3.0) 

27 

( 6 . 0 ) 


222-236 

(50.0-53.0) 

276 

(62.0) 

311-316 

(70.0-71.0) 

356 

(80.0) 

400 

(90.0) 

449 

( 100 . 0 ) 

489 

( 110 . 0 ) 

525 

(118.0) 

569 

(128.0) 


13-18 
(3. 0-4.0) 

18 

(4.0) 

13 

(3.0) 

13 

(3.0) 

18 

(4.0) 

18 

(4.0) 

22 

(5.0) 

18 

(4.0) 

27 

( 6 . 0 ) 


14 145 


5 000 


5 229 


3 161 


1 798 


1 402 


5.685 


5.834 


6.152 


6.459 


6.779 


7.165 


7.369 


7.564 


7.724 


2.238 


2.297 


2.422 


2.543 


2.669 


2.821 


2.901 


2.978 


3.041 


Total 


31 970 


0O 

Ui 

















Nominal 

Horizontal load, 


P„ , kN 
Max 

P 

H > kN (kip) 


(kip) 

High max 

Low max 

Min 


LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-5 


Vertical load* Increment 


(44.4) 

205 

(46.0) 

225 

(50.6) 

229 

(51.4) 

254 

(57.0) 


P y , kN (kip) 


High max 


(33.6-34.9) 

177 

(39.8) 

193 

(43.4) 

215 

(48.4) 

243 

(54.6) 

240 

(54.6) 

265 “ 

(59.5) 


Increment 
of cycles, 
4N 


39-41 
(8. 8-9.3) 

51-54 

(11.4-12.2) 

54-55 

(12.2-12.4) 

69 

(15.6) 

73 

(16.4) 

84 

(18.8) 

98 

( 22 . 0 ) 

96 ^ 

( 21 . 6 ) 

113 

(25.5) 


3-5 

25 264 

(0.6-1. 2) 


6-8 

6 425 

(1.4-1. 8) 


3-8 

5 462 

(0.6-1. 8) 



Crack length, 2a 


Horizontal 

crack 

cm (in.) 


2a. = 5.283 
= (2.080) 

5.550 

(2.185) 

5.923 

(2.332) 

6.500 

(2.559) 

7.059 

(2.779) 

7.529 

(2.964) 

8.082 

(3.182) 

8.506 

(3.349) 

8.948 

(3.523) 

9.373 

(3.690) 


Vertical 

crack 

cm (in.) 


5. '*74 
(2.155) 

5.474 

(2.155) 

5.608 

(2.208) 

5.756 

(2.266) 

5.946 

(2.341) 

6.093 

(2.399) 

6.185 

(2.435) 

6.210 

(2.445) 

6.231 

(2.453) 


43 105 
















LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-6 


Nominal 

P M » kN 

Max - 

Horizontal load, 

T-S 1 IT /I « \ 

' Vertical load, 

P , kN (kip) 

V 

Increment 
of cycles, 
AN 

Crack length, 2a 

kin 

<Kip; 

cm 

in* 

(kip) 

Max 

Min 

Max 

Min 





/ 




2a. - 4.928 

l 

2a. = 1.940 

i 

116 

(26) 

120 

(27.0) 

17 

(3.8) 

121 
('27. 2) 

17 

(3.8) 


5.197 

2.046 

138 

(31) 

141 

(31.6) 

17 

(3.8) 

137 

(30.8) 

15 

(3.4) 

15 000 

6.139 

2.417 

160 

(36) 

156 

(35.0) 

18 

(4.0) 

157 

(35.2) 

16 

(3.6) 

5 013 

6.614 

2.604 


155 

(34.8) 

18 

(4.0) 

157 

(35.2) 

13 

(3.0) 

2 917 

6.914 

2.722 


157 

(35.2) 

18 

(4.0) 

157 

(35.4) 

16 

(3.6) 

5 000 

7.209 

2.838- 

178 

(40) 

179 

(40.2) 

19 

(4.2) 

180 

(40.4) 

19 

(4.2) 

1 510 

7.523 

2.962 


180 

(40.4) 

20 

(4.6) 

181 

(40.6) 

21 

(4.8) 


7.851 

3.091 


180 

(40.4) 

20 

(4.6) 

181 

(40.6) 

21 

(4.8) 


8.197 

3.227 

200 

(45) 

199 

(44.8) 

24 

(5.4) 

198 

(44.6) 

20 

(4.4) 

511 

8.395 

3.305 


200 

(45.0) 

24 

(5.4) 

198 

(44.6) 

20 

(4.6) 


8.65_4_ 

3.407 


199- 

(44.8) 

22 

(5.0) 

199 

(44.8) 

19 

(4.2) 

1 000 

9.352 

3.682 

222 

(50) 

222 

(49.8) 

25 

(5.6) 

221 

(49.6) 

26 

(5.8) 

509 

9.873 

3.887 


222 

(49.8) 

25 

(5.6) 

218 

(49.0) 

23 

(5.2) 


10.36 

4.080 


222 

(50.0) 

25 

(5.6) 

218 

(49.0) 

24 

(5.4) 


10.91 

4.297 

245 

(55) 

244 

(54.8) 

25 

(5.6) 

237 

(53.2) 

24 

(5.4) 

210 

11.37 

4.478 


242 

(54.4) 

\ 

25 

(5.6) • 

244 

(54.8) 

28 

(6.4) 


11.82 

j 

4.652 


243 

(54.6) 

24 

(5.4) 

242 

(54.4) 

28 

(6.2) 


12.27 

4.832 





Total 

55 005 




187 














LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-7 


Nominal 

Horizontal load. 

Vertical load, 

Increment 

Crack length, 

P„ > kH 
Max 

P R , kN (kip) 

P v , kN (kip) 

of 

cycles. 

2a 


(kip) 






AN* 



Max 

Min 

Max 

Min 


cm 

in* 








2a. = 5.166 
i 

2.034 

289 

.282 

24 

287 

31 

4 982 

5.527 

2.176 

(65) 

(63.5) 

(5.5) 

(64.5) 

(7.0) 

5 

000 

5.878 

2.314 






5 

ooo 


2.449 

311 

311 

29 

322 

31 

3 

467 

6.642 

2.615 

(70) 

(70.0) 

(6.5) 

(72.5) 

(7.0) 

3 

500 

7.046 . 

• 2.774 






3 

500 

7.523 

2.962 

343 

347 

31 

345 

38 

2 

331 

7.925 

3.120 

(77) 

(78.0) 

(7.0) 

(84.5) 

(8.5) 

2 

000 

8.351 

3.288 






2 

000 

8.745 

3.443 

378 

376 

31 

376 

38 

1 

213 

9.103 

3.584 

(85) 

(84.5) 

(7.0) 

(84.5) 

(8.5) 

1 

200 

9.535 

3.754 






1 

200 

9.837 

3.873 

409 

409 

36 

409 

42 


755 

10.23 


(92) 

(92.0) 

(8.0) 

(92.0) 

(9.5) 


750 

10.68 

4.205 







750 

11.14 

4.387 







750 

11.68 

4.600 

445 

440 

44 

445 

44 


306 

12.17 

4.791 

(100) 

(99.0) 

(10.0) 

(100.0) 

(10.0) 


300 

12.54 

4.937 







300 

12.92 

5.088 

476 

472 

49 

476 

49 


155 

13.36 

5.259 

(107) 

(106.0) 

(11.0) 

(107.0) 

(11.0) 


150 

13.84 

5.449 







150 

14.31 

5.633 

512 

507 

58 

512 

53 


80 

14.75 

5.805 

(115) 

(114.0) 

(13.0) 

(115.0) 

(12.0) 


75 

15.31 




i 




75 

15.91 

6.264 







75 

16.61 

6.538 







75 

17.42 

6.860 







75 

18.33 

7.218 







75 

19.70 

7.757 







75 

21.54 

8.479 







75 

24.98 

9.833 

Fracture 

503 


516 



16 

28. 88 a 



(113.0) 


(116.0) . 







Total 

40 

455 

11.37 a 


3 . 

Extrapolated from compliance data. 


188 
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LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-8 


Nominal 

, kN 
Max’ 

(kip) 

Horizontal load, 
P H , kN (kip) 

Vertical load, 

Increment 
of cycles,, 
AN 


Crack length, 

2a 

V 


Horizontal 

Vertical 

High 

Low 

Min 

High 

Low 

Min 

crack 

crack 


max 

max 


max 

max 



cm 

in. 

cm m. 









2a. = 6.375 
. i 

2a. = 2.510 


289 

(65) 

287 

(64.5) 

in 

(25.0) 

16 

(3.5) 

276 

(62.0) 

107 

(24.0) 

18 

(4.0) 

346 

6.584 

2.592 

Disregard 

289 

(65) 

285 

(64.0) 

107 

(24.0) 

11 

(2.5) 

276 

(62.0) 

105 

(23.5) 

13 

(3.0) 

6 525 . 

6.965 

2.742 


289 

'(65) 

285 

(64.0) 

109 

(24.5) 

13 

(3.0) 

276 

(62.0) 

107 

(24.0) 

16 

(3.5) 

6 000 

7.602 

2.993 


320 

(72) 

311 

(70.0) 

122 

(27.5) 

16 

(3.5) 

314 

(70.5) 

129 

(29.0) 

24 

(5.5) 

3 006 

7.950 

3.130 


351 

(79) 

349 

(78.5) 

140 

(31.5) 

16 

(3.5) 

343 

(77.0) 

140 

(31.5) 

22 

(5.0) 

1 500 

8,235 

3.242 


383 

(86) 

378 

(85.0) 

156 

(35.0) 

20 

(4.5) 

365 

(82.0) 

151 

(34.0) 

22 

(5.0) 

1 500 

8.580 

3.378 


383 

(86) 

378 

(85.0) 

156 

(35.0) 

m 

365 

(82.0) 


22 

(5.0) 

1 500 

9.055 

3.565 


414 

(93) 

413 

(93.0) 

169 

(38.0) 

36 

(8.0) 

427 

(96.0) 

182 

(41.0) 

36 

(8.0) 

750 

9.444 

3.718 


445 

(100) 

423 

(95.0) 

187 

(42.0) 

36 

(8.0) 

418 

(94.0) 

182 

(41.0) 

31 

(7.0) 

750 

9.779 

3.850 


476 

(107) 

458 

(103.0) 

185 

(41.5) 

44 

(10.0) 


189 

(42.5) 

33 

(7.5) 

514 

10.33 

4.065 


476 

(107) 

458 

(103.0) 

185 

(41.5) 

44 

(10.0) 


189 

(42.5) 

33 

(7.5) 

500 

10,95 

4.309 


507 

(114) 

485 

(109.0) 

205 

(46.0> 

49 

(11.0) 

478 

(107.5) 

205 

(46.0) 

33 

(7.5) 

350 

11.83 

4.659 


507 

(114) 

480 

(108.0) 


53 

(12.0) 


205 

(46.0) 

31 

(7.0) 

169 

12.42 

4.890 


538 

(121) 

514 

(115.5) 

222 

(50.0) 

47 

(10.5) 


222 

(50.0) 

40 

(9.0) 


13.17 

5.185 


538 

(121) 

512 

(115.0) 

218 

(49.0) 

47 

(10.5) 


220 

(49.5) 

38 

(8.5) 


13.83 

5.455 


538 

(121) 







63 

14.73 

q 7QR 








j • / yo 








Total 

23 673 



• 






LOADING SUHEDULL AND GKAGK GKUWIH £ UK GIGL1G TEST OK SN-y 



Extrapolated from crack growth rate. 























LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-27 


Nominal 
P * kN 

Horizontal load, 

P . IcN f'ln'n'i 

Vertical load, 

Increment 
of cycles, 
AN 

Crack length, 2a 

Max’ 

(kip) 

H 

N IT J 



cm 

in. 

Max 

Min 

Max 

Min 









2a. = 7.701 

i 

2a. = 3.032 

i 

267 

(60) 

269 

(60.5) 

20 

(4.5) 

265 

(59.5) 

18 

(4.0) 

884 

7.904 

3.112 


267 

(60.0) 

20 

(4.5) 

258 

(58.0) 

20 

(4.5) 

5 481 

8.197 

3.227 


267 

(60.0) 

18 

(4.0) 

269 

(60.5) 

18 

(4.0) 

4 970 

9.195 

3.620 

267 

(60) 

269 

(60.5) 

20 

. (4.5) 

274 

(61.5) 

24 

(5.5) 

1 732 

9.614 

3.785 


267 

(60.0) 

18 

(4.0) 

265 

(59.5) 

20 

(4.5) 

1 597 

9.962 

3.922 


267 

(60.0) 

18 

(4.0) 

267 

(60.0) 

18 

(4.0) 

1.719 

10.34 

4.081 

267 

(60) 

267 

(60.0) 

18 

(4.0) 

267 

(60.0) 

20 

(4.5) 

1 442 

10.78 

4.245 


271 

(61.0) 

18 

(4.0) 

267 

(60.0) 

18 

(4i0) 

1 788 

11.30 

4.447 


271 

(61.0) 

18 

(4.0) 

262 

(59.0) 

18 

(4.0) 

1 307 

11.66 

4.590 

311 

(70) 

316 

(71.0) 

.22 

(5.0) 

309 

(69.5) 

20 

(4.5) 

786 

12.11 

4.768 


318 

(71.5) 

24 

(5.5) 

311 

(70.0) 

22 

(5.0) 

895 

12.66 

4.983 


318 

(71.5) 

27 

(6.0) 

309 

(69.5) 

20 

(4.5) 

899 

13.20 

5.198 

356 

(80) 

360 

(81.0) 

22 

(5.0) 

356 

(80.0) 

24 

(5.5) 

408 

13.90 

5.471 


363 

(81.5) 

22 

(5.0) 

354 

(79.5) 

22 

(5.0) 

354 

14.40 

5.671 


360 

(81.0) 

22 

(5.0) 

356 

(80.0) 

27 

(6.0) 

392 

14.95 

5.884 

400 

(90) 

405 

(91.0) 

29 

(6.5) 

387 

(87.0) 

27 

(6.0) 

160 

15.44 

6.079 


407 

(91.5) 

27 

(6.0) 

387 

(87.0) 

27 

(6.0) 

150 

15.98 

6.290 


407 

(91.5) 

29 

(6.5) 

387 

(87.0) 

24 

(5.5) 

150 

16.40 

6.458 

445 

(100) 

431 

(97.0) 

44 

(10.0) 

440 

(99.0) 

31 

(7.0) 

72 

16.84 

6.629 


436 

(98.0) 

44 

(10.0) 

436 

(98.0) 

31 

(7.0) 

60 

17 . 16 

6.757 


436 

(98.0) 

44 

(10.0) 

436 

(98.0) 

27 

(6.0) 

60 

17 . 42 

i 

6.860 

445 

(100) 

436 

(98.0) 

44 

(10.0) 

436 

(98.0) 

31 

(7.0) 

120 

18.09 

7.116 


436 

(98.0) 

44 

(10.0) 

436 

(98.0) 

36 

(8.0) 

120 

18.81 

7.404 


436 

(98.0) 

44 

(10.0) 

436 

(98.0) 

31 

(7.0) 

120 

19.56 

7.700 




LOADING SCHEDULE AND CRACK GROWTH FOR CYCLIC TEST OF SN-27 (concl) 


Nominal 

Max 

(kip) 

Horizontal load, 
Pg) kN (kip) 

Vertical load, 
P v , kN (kip) 

Increment 
of cycles , 
AN 

Crack length, 2a 

cm 

in* 

Max 

Min 

Max 

'Min 

445 

436 

Kami 

436 

31 

120 

20.53 

8.083 


(98.0) 


(98.0) 

(7.0) 





431 

1 

436 

27 

120 

21.43 

8.436 


(97.0) 

1 

(98.0) 

(6.0) 





436 

44 

436 

27 

120 

22.57 

8.886 


(98.0) 

(10.0) 

(98.0) 

(6.0) 




445 

436 

44 

436 

27 

. 120 

23.93 

9.421 

(100) 

(98.0) 

(10.0) 

(98.0) 

(6.0) 





436 

AQ 

436 

27 

60 




(98.0) 


(98.0) 

(6.0) 




1 



Total 

26 206 




192 



















APPENDIX E 


PHOTOGRAPHS OF SPECIMENS AFTER TESTING 


Note: This appendix contains photographs of all the specimens 

after testing, in numerical order. 



Figure 1. - SN-1 after Fracture, Console Side. 



Figure 2. - SN-1 after Fracture, Opposite Side. 
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Figure 3. - Closeup of Console Side, SN-1 . 


« 



Figure 4. - Closeup of Opposite Side, SN-1. 
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Figure 5. - Overall View of SN-2 after Testing, Console Side. 



Figure 6. - Overall View of SN-2 after Testing, Opposite Side. 
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Figure 7. - Closeup View of SN-2 after Fracture, Console Side. 
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Figure 8. - Closeup View of SN-2 after Fracture, Opposite Side 
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Figure 10. - Central Area of SN-3 after Testing, Opposite Side. 
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Figure 11. - Closeup of SN-3 after Testing, Console Side. 



Figure 12. - Closeup of SN-3 after Testing, 
Opposite Side. 
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Figure 14. - Overall View of SN-4 after Testing, Opposite Side. 
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Figure 15. - Closeup of SN-4 after Testing, Console Side, 



Figure 16. - Closeup of SN-4 after Testing, Opposite Side. 
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Figure 17. - Overall View of SN-5 after Testing, 
Console Side . 



Figure 18. - Overall View of SN-5 after Testing 
Opposite Side. 


» 
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Figure 19. - Closeup of SN-5 after Cutting, Console Side. 



Figure 20. - Closeup of SN-5 after Cutting, Opposite Side. 
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Figure 21. - Overall View of SN-6 after Testing, 
Console Side . 
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Figure 23. - Closeup of SN-6 after Testing, 
Console Side . 



Figure 24. - Closeup of SN-6 after Testing, Opposite Side. 
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Figure 26. 


- Overall View of SN-7 after Testing, 
Console Side . 


- Overall View of SN-7 after Testing, 
Opposite Side . 


Figure 25. 
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206 



Figure 27. - Closeup of SN-7 after Testing, 
Console Side . 



Figure 28. - Closeup of SN-7 after Testing, 
Opposite Side . 
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Figure 29. - Overall View of SN-8 after Testing. 
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Figure 30. - Central Portion of SN-8 after Testing. 
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Figure 31. - Overall View of SN-9 after Testing, Console Side. 



Figure 32. - Overall View of SN-9 after Testing, Opposite Side. 
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Figure 33. - Central Portion of SN-9 after Testing, 
Console Side. 



Opposite Side. 
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Figure 35. - SN-11 after Test. 



Figure 36. - Central Portion of SN-11 after Test. 
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Figure 38. - SN-12 after Testing. 
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Figure 39. - Central Portion of SN-12 
after Testing. 



Figure 40. - SN-14 after Fracture Test. 
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Figure 41. - Central Portion of SN-14 after Test, Console Side. 



Figure 42. - Central Portion of SN-14 after Test, Opposite Side. 
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Figure 43. - Overall View of SN-15 after Testing. 



Figure 44. - Central Portion of SN-15 after Testing. 
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Figure 45. - Overall View of SN-16 after Testing. 
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Figure 47. - SN-17 after Testing. 



Figure 48. - Central Portion of SN-17 after Testing. 
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Figure 51. - SN-18 after Test, Central Region. 



Figure 52. - Fracture Surface of SN-18. 
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Figure 53. - SN-23 after Test. 



Figure 54. - SN-23 after Test, Central Region. 
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Figure 55. - SN-24 after Test . 



Figure 56. - SN-24 after Test, Central Region. 
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Figure 59. - SN-25 after Test, Central Region. 



Figure 60. - Fracture Surface of SN-25. 
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Figure 62. - SN-26 after Test, Central Region. 
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Figure 63. - SN-27 after Testing. 



Figure 64. - Central Portion of SN-27 after Testing. 
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APPENDIX F SYMBOLS 


A 


a 


a 


cr 


a, 

1 


a, . 


11 



H 

a. 

1 


I 


a . 

l 


ia VB 

Aa VT 


Area of specimen appropriate to applied stress in the 
expressions for stress intensity* 

M 

One-half actual crack length* 

Critical value of a* 

One-half initial crack length, after precrackirig. 
Coefficients in Sih theory* 


a^ [2E/ (1+v) (1-2v) ] 

One-half the arc length of curved flaws after cyclic, but 

before static, testing. Used for Hols ton analysis. 

\ 

The projection of a^ perpendicular to final fracture di- 
rection. Used for' Iida analysis. Also, one half the 
projection of branched cracks after cyclic but before 
static testing. 


Increment of branched crack growth at the bottom in the 
horizontal direction. 


Increment of branched crack 
horizontal direction. 

Increment of branched crack 
vertical direction . 

i 

Increment of branched crack 
vertical direction. 


growth at the top in the 
growth at the bottom in the 
growth at the top in the 


b One-half plate width — Iida analysis. 

b^ Six times the equivalent linearized crack slant length 

for the Iida analysis. Used to establish relation between 
crack length and equivalent plate width of Maltese cross 
specimens. 

C Compression. Also generalized coefficient in stress in- 

tensity expression. 

Generalized coefficient in crack growth equation. 
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COD 

CSD 

E 

G 

G 

cr 

K. 

1 

k. 

i 

K t 


K 


II 


K 


Ic 


K 


He- 


Iu 


K. 

K, 


IIu 


II 


M. 


SI 


M, 


SII 


N 

AN 


n 


P 

P 


Crack opening displacement. 

Crack shear displacement. 

Young’s modulus. 

Crack extension force. 

Critical value of G.' 

Stress intensity in the ith mode. 



Mode I stress intensity. 

Mode II stress intensity. 

Critical value of k^. 

Critical value of K^. 

.Mode I stress intensity at fracture. 

Mode II stress intensity at fracture. 

Mode I stress intensity coefficient" in Holston analysis. 

Mode II stress intensity coefficient in Holston analysis. 
Measure of plastic zone size.. . . . . 

Mode I stress intensity coefficient in lida analysis. 

Mode II stress intensity coefficient in lida analysis. 
Number of cycles, of loading. 

Increment of. cycles per load block.. . 

Experimental value of at any time. Also curve 

fitting exponent in crack growth equation.- 

Load., " • ' " •• * 

Load perpendicular to predominant crack direction or 
branch, .lida analysis-. ' •: 


228 



^3 S CD CD ~ r< P <? P 


H 


V 


R 

S 

T 

t 

U 


a 

a 


a 



a 


c 



Load in A axis direction. 

Load in B axis direction. 

Horizontal load. 

Vertical load. 

Radius of arc crack. Also crack extension resistance. 
Intensity of strain energy density, Sih theory. 

Tension. Also period of cyclic load. 

Thickness. Also time lag of horizontal ► load axis. 

Strain energy density, Sih theory. 

Crack shearing displacement normalized to infinite plate 
value . 

Crack opening displacement normalized to infinite plate 
value . 

One-half interior angle of sector defined by arc crack. 
Arctan 

Angle in polar coordinates defining the predicted fracture 
direction. 

a/b= normalized projected crack length. 

Poisson's ratio. 

v for plane strain; v/(l+v) for plane stress. 

Normal stress. 

Sress parallel to load direction. 

Stress perpendicular to load direction. 

Stress at center of unflawed plate under equal biaxial 
loads . 

Critical value of o. 

Horizontal stress. 
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Assumed stress levels for construction of R curve, 

3. 1 5 2, ♦ ♦ « 

Normal stress. 

Vertical stress. 

Yield strength. 

Tangential stress in polar coordinate system. 

Shear stress. 

Arctan 

Arctan 31^ - 
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